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SUMMARY 


This  is  a research  program  that  seeks  to  incorporate  EEG  "evoked  responses" 
to  complement  overt  behavior  in  specific  man-machine  communication  schemes 
involving  various  aspects  of  decision  making.  Behind  the  experimental  strategy 
is  the  view  that  the  EEG  in  evoked  responses  is  made  up  of  a complex  of  wave- 
lets that  reflect  individual  and  sequential  events  in  the  brain  and  in  particular  the 
cerebral  cortex.  T$ie  rules  for  appearance  of  these  elementary  wavelets  would 
be  the  "syntactic"  constraints  upon  a neuroelectric  language  with  the  wavelets 
acting  as  individual  signatures  for  various  aspects  of  the  neural  transactions 
that  follow  presentation  of  a stimulus.  Once  identified,  the  reliability  and 
stability  of  those  signatures  can  be  reinforced  naturally  in  the  man-machine 
communication  scheme,  in  a way  amounting  to  operant  conditioning.  Thus  the 
approach  is  centered  on  sequential  events  of  short  duration  in  bio-electric 
potentials  and  the  relation  between  these  sequences,  behaviors  and  brain  states. 
Those  could  in  fact  lead  to  a new  model  for  the  phenomena  of  evoked  responses 
in  the  EEG. 

Current  experiments  are  aimed  at  visual  evoked  responses  using  colored  and 
patterned  visual  stimuli,  void  at  first  of  any  cognitive  or  emotional  content. 
Color  flashes  and  patterns  are  choice  stimuli  that  would  provide  potential  sup- 
port for  non-verbal  symbols  in  a man-machine  communication  language.  In 
the  color  experiments,  both  target  and  intensities  in  each  color  are  varied 
separately.  From  the  ensemble  of  these  experiments,  tentative  candidate  EEG 
codes  emerged  that  appear  to  reflect  three  retinal  processes  (red-green-blue). 
Best  results  have  been  obtained  with  a yellow  background  apparently  because  of 
the  red-green  depression  caused  by  retinal  adaption  to  yellow  light.  This 
allows  the  blue  component,  which  had  been  elusive  before,  to  appear  much  more 
clearly.  These  are  remarkable  results  since,  while  there  is  good  evidence  for 
a trichromatic  absorbing  structure  in  the  fovea,  there  was  no  indication  that 
these  processes  would  be  reflected  into  scalp  potentials. 


1 


Topological  Dimensionality  studies  of  the  same  data  showed  that  red  and 
green  components,  although  close  in  timing,  did  behave  in  a clearly  distinct 
manner.  In  addition  the  response  of  the  dark  adapted  eye,  a limit  case  in 
background  levels,  was  found  clearly  detached  from  the  cluster  in  all  cases. 
Simultaneously  multichannel  single  epochs  from  the  same  experiments,  have 
been  analyzed  off-line  using  stepwise  discriminant  analysis  with  and  without 
the  detour  of  an  orthogonal  transformation.  The  results  have  been  startling; 
the  stimulus  sets  have  proven  extremely  effective  in  providing  a quality  of 
discrimination  that  was  unheard  of  with  single  EEG  epochs,  using  only  ten 
amplitude  samples.  Various  representations  were  attempted,  among  which 
the  replacement  of  the  single  evoked  response  by  the  point  processes  formed 
by  the  instants  at  which  a single  epoch  was  displaying  a peak  in  the  positive 
or  negative  direction.  Cumulative  histograms  of  those  micro-events  were 
found  to  be  as  stimulus -specific  as  are  the  averages,  and  it  appears  that 
some  of  the  observed  peaks  correspond  to  the  firing  of  specific  nuclei.  This 
raises  great  hope  that  a combination  of  short  event  detection  with  amplitudes 
at  selected  times  will  provide  encoding  of  the  evoked  response  in  a particularly 
economical  manner  in  man-machine  communication.  Already  at  this  time, 
even  without  the  benefit  of  amplitude  information  the  peak  "message"  has 
equivalent  clarity  or  discriminability,  when  compared  to  the  peak  histogram, 
than  the  single  epoch  with  respect  to  the  average. 

The  experiments  have  raised  continuous  challenges  to  the  computer  system 
and  require  additional  developments.  The  central  software  package  has  been 
successfully  implemented  on  the  IBM  360/91  at  UCLA  and  will  be  accessible 
through  the  ARPANET.  It  consists  of  a sequence  of  "code  cracking*'  programs 
operating  on  epoch  oriented  data  that  will  be  ultimately  used  interactively. 
Several  hardware  improvements  were  made  on  the  indigenous  laboratory  com- 
puters and  the  high-speed  link  that  connects  the  laboratory  to  the  360/91. 

A tape  drive  and  line  printer  controller  were  added  to  the  SDS  920.  Long 
vector  hardware  was  successfully  retrofitted  on  the  IMLAC  PDS-1  to  allow 
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the  creation  of  pattern  stimuli  under  computer  control. 


Following  recent  directives,  the  whole  system  has  been  re-evaluated  in  view 
of  its  potential  role  in  the  program  of  data  sharing  and  interconnection  between 
Biocybernetic  Laboratories.  Special  attention  has  been  given  to  the  articulation 
of  the  BCI  software  package  with  software  resources  at  MIT-MULTICS.  Finally 
a remote  hardwired  terminal  connection  between  CCBS  and  BCI  has  been  imple- 
mented and  is  now  in  operation. 
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HUMAN  EVOKED  RESPONSE  EXPERIMENTS 


In  1929  Berger  demonstrated  the  possibility  of  recording  brain 
waves  from  the  intact  skull.  Since  then,  an  enormous  amount 
of  brain  wave  data  covering  a variety  of  conditions  has  been 
accumulated  by  neurophysiologists,  and  in  recent  years,  com- 
puters have  been  used  extensively  for  analysis. 

Overall  characteristics  of  these  fluctuations  of  electrical  poten- 
tial can  be  somewhat  predicted  in  relation  to  the  electrode  site, 
the  mental  state  of  the  subject,  and  the  presence  and  type  of  sen- 
sory stimulation.  Some  of  those  characteristics  are  readily 
identified  by  eye.  Well-known  examples  are  recognition  of  alpha 
activity  and  the  phenomenon  of  alpha  blocking,  sleep  and  barbi- 
tuate  spindles,  and  the  3 -per -second  spike  and  wave  complex  of 
petit  mal  epilepsy.  More  subtle  information  in  the  EEG  signals, 
however,  requires  computer  analysis. 


This  "spontaneous"  or  "on-going"  electrical  activity  is  somewhat 
rhythmic  in  nature.  The  analysis  of  these  rhythms  has  retained 
much  of  the  early  attention  paid  to  brain  waves  in  general.  Yet 
it  is  increasingly  evident  that  rhythmic  activities  contain  little 
information  in  themselves  and  that  their  function  is  probably 
similar  to  that  of  a carrier.  The  general  picture  is  that  idle 
nervous  tissue  will  exhibit  spontaneous  oscillation  or  rhythm  while 
activity  or  commitment  of  the  same  tissue  to  an  active  function 
will  be  denoted  by  desynchronized  random-like  oscillations. 
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Beyond  gross  differentiation  of  brain  states,  it  seems  that 
information  coding  in  the  EEG  wave  should  be  sought  in  the 
specific  waveforms  generated  in  time.  By  contrast  with  the 
wide  spread  character  of  oscillations,  the  wave  shapes  obtained 
with  "evoked  potential"  are  localized  and  correspond  well  with 
underlying  post-synaptic  potentials.  Evoked  potentials 
(evoked  responses)  are  generally  obtained  with  adight  flash, 
brief  sound,  or  touch  of  the  skin  which  gene  rate  ^in  the  cor- 
responding sensory  cortex  (visual,  auditory  or  somesthetic) , 
a localized  electrical  response  in  the  form  of  an  aperiodic  wave 
form  lasting  up  to  half  a second  and  superimposed  on  the  on- 
going background  activity.  In  general,  repeated  stimuli  and 
averaging  of  the  waveforms  have  been  used  to  reveal  the 
"evoked"  response  by  cancelling  the  background  "noise". 


Indeed  early  investigations  in  1965  by  Fox  and  O'Brien  resulted 
in  the  demonstration  that  stimulus-evoked  potentials  were 
directly  related  to  the  probability  of  firing  of  any  particular 
cell  in  the  area  of  population  included  in  the  brain  wave  recording. 
Specifically,  the  studies  showed  that  in  these  cortical  cells, 
the  probability  of  firing  could  be  described  almost  precisely  by 
the  shape  of  the  slow  wave  associated  with  that  same  electrode. 
Thus,  the  relationship  between  single -cell  spikes  and  evoked 
potential  was  dramatically  clarified. 


Fox  and  Norman  demonstrated  that  the  spike -wave  relationship 
would  also  hold  for  spontaneous  activity  of  the  cortex.  The  next 
step  then  was  to  determine  if  these  potentials  might  represent 
behavior,  based  on  these  earlier  studies  which  showed  that  the 
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moment -to -moment  changes  in  amplitude  and  polarity  of  the 
slow  activity  of  the  brain  represented  moment-to -moment 
changes  in  probability  of  firing  of  single  cells. 

To  achieve  this  took  a novel  approach.  Fox  and  Rudell  for  the 
first  time  used  operant  reinforcement  techniques  to  increase  or 
decrease  the  probability  of  occurrence  of  some  component  of 
the  visual  evoked  potential.  This  they  claimed  was  tantamount 
to  asking  the  animal  whether,  under  reinforcement  control,  it 
could  increase  or  decrease  the  probability  of  some  aspect  of  his 
brain  wave.  In  turn  this  was  equivalent  to  asking  whether  or 
not  the  animal  had  any  behavior  available  for  reinforcement, 
which  utilizes,  encodes  or  is  represented  by  the  particular  aspect 
of  brain  activity  chosen. 

To  pursue  the  study  of  this  brain-behavior  relationship  to  time- 
dependent  behaviors,  Rosenfeld  and  Fox  trained  unrestrained 
cats  with  implanted  electrodes  to  make  discrete  paw  movements. 

At  the  same  time  the  mixed  sensory-motor  potential  was  recorded 
in  the  cortex.  These  data,  treated  correlatively , showed  a relation 
ship  of  early  waves  and  early  portions  of  each  movement  and  of 
later  waves  and  later  portions  of  each  movement,  but  relative 
independence  of  early  and  late  measurements. 

In  the  second  part  of  these  experiments  a particular  aspect  of  the 
movement-evoked  potential  in  the  cortex  was  chosen  and  by 
reinforcement  control,  the  animals  were  trained  to  increase  the 
probability  or  amplitude  of  this  particular  component  of  the  brain 
wave.  Exactly  in  accordance  with  prediction,  the  animals,  when 


this  training  was  accomplished,  altered  some  aspect  of  their 
limb  movements.  To  the  discrete  changes  of  the  cortical 
evoked  potential  associated  with  movement, correspondea  some 
discrete  and  finely  detailed  changes  in  some  aspects  of  the  move- 
ment of  the  limb,  showing  that  spontaneous  activity  of  the  brain 
from  moment-'co -moment  encodes  function  and  the  sequential 
probabilities  of  single -cell  firing  were  functionally  represented 
by  the  sequential  and  momentary  changes  in  brain  waves.  More 
important,  moment-to-momsnt  changes  in  behavior  were  dis- 
cretely reflected  in  moment-to-moment  changes  in  sequential 
brain  waves . 

These  findings  in  animal  studies  are  very  significant  insofar  as 
they  point  out  the  functional  meaning  of  the  "gross"  EEG  potentials 
generated  by  populations  rather  than  single  neurons.  In  other 
words,  at  leaBt  in  some  conditions,  the  neural  "roar"  does  not 
cancel  out  into  hopeless  noise.  With  humans  and  without  im- 
planted electrodes , the  difficulties  are  certainly  compounded. 

But  in  compensation  more  avenues  are  open  for  operant  condi- 
tioning, a technique  that  may  well  be  the  key  to  success  in  identi- 
fying the  EEG  "codes". 


A starting  point  in  this  study  was  the  recognition  that  different 
stimuli  with  specific  set  of  features  evoke  distinguishable  electri- 
cal "signatures"  on  the  scalp.  For  instance,  the  response  to  the 
brief  hashing  of  a figure  made  of  vertical  lines  will  yield  a wave- 
form markedly  different  from  that  obtained  from  a set  of  circles. 

In  fact,  the  presence  in  the  evoked  waveforms  of  clear  correlates  of 
the  modalities  of  sensory  stimulation  has  been  abundantly  demonstrated. 
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Of  particular  interest  are  studies  dealing  with  visual  stimuli 
{White  and  Eason,  1966;  Karter  and  White,  1968,  1969;  Clynes 
and  Kohn,  1967;  Rietveld,  et  al,  1967;  Spehlmann,  1965; 
Spekreise,  1966).  The  evoked  electrical  signature  on  the 
cortex  does  however  include  more  than  a mirror  conversion 
of  the  stimulus  content.  In  fact  only  the  early  part  of  the  response 
appears  to  be  directly  stimulus  bound,  while  the  "late"  com- 
ponents appear  to  relate  to  more  complex  brain  functions  re- 
lated to  the  stimulus  such  ais  its  perception  and  meaning. 

The  general  approach  was  defined  as  follows  in  the  original 
proposal:  First,  stimulus  sets  are  selected  that  span  some 
"dimensions"  expected  to  elicit  neural  correlates,  that  is, 
changes  in  evoked  response.  Next,  wave  components  in 
the  neural  response  are  identified  that  appear  to  correlate 
with  each  stimulus  dimension,  thus  providing  a first  identi- 
fication of  the  "candidate"  code,  along  with  its  sensitivity 
and  stability.  In  this  phase  the  model-fTee  topological  dimen- 
sionality analysis  is  applied  to  the  data  to  reveal  its  parametric 
dimensionality  in  relation  to  each  stimulus  dimension.  Out- 
come "candidate"  waveforms  or  parameters  are  then  isolated 
and  compared  from  subject  to  subject.  For  each  subject  the 
reliability  and  stability  of  the  emerging  dimensions  must  be  eval- 
uated in  repeated  experiments.  As  the  initial  bioelectric  con- 
figuration emerges  from  dimensionality  analysis,  an  interactive 
closed  loop  mode  is  entered.  Operant  conditioning  techniques 
are  applied  to  individualize  and  sharpen  these  waveforms. 
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• 

Current  experiments  are  aimed  at  the  conditioning  of  selected  attributes 

• of  visual  evoked  responses  and  are  using  colored  and  patterned  visual 
stimuli  (void  at  first  of  any  cognitive  or  emotional  content).  Color 
flashes  and  patterns  are  choice  stimuli  that  would  provide  potential  support 
for  non-verbal  symbols  In  a man-machine  communication  language. 

• Eventually  the  cognitive  or  emotional  content  of  the  stimulus  was  to  be 
considered  (e.  g.  whether  the  subject  has  directed  his  attention  towards 
a specific  feature  in  the  complex  stimulus,  or  whether  the  occurrence 

• of  the  stimulus  represents  a gain  or  a loss  In  a game  playing  situatlon- 
or  whether  the  stimulus  Is  varied  in  a dimension  varying  from  pleasant 

I to  aversive).  These  extensions  are  to  be  attempted  during  the  second 

and  third  year. 

$ 

A brief  review  of  the  basic  laws  of  color  vision  will  be  given  next  as 
I an  introduction  to  the  experiments  on  evoked  response  to  color. 

k © 
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TRICHROMATICITY 


The  essential  trichromaticity  of  color  vision  has  been  demonstrated  by 
tile  results  of  oolor  mixing  experiments  using  psychophysical  technique 
for  making  color  judgment,  i.e.  using  the  eye  as  a null  instrument. 

Any  spectral  color  can  be  matched  by  the  proper  mixture  of  any  three 
colors,  as  long  as  none  of  those  is  a combination  of  the  other  two. 

A color  stimulus  can  be  characterized  by  its  C.I. E.  (Commission  Inter- 
national de  l'Eclairage)  tristimulus  values  (X,  Y,  Z)  i.e.  by  the  relative 
amount  of  three  standard  C.  I.  E.  priorities  that  have  been  required  by 
the  "standard  observer"  to  match  the  stimulus.  Tristimulus  values  and 
luminosity  functions  (visual  efficiency)  in  function  of  wavelength  (ire.  for 
pure  colors)  are  given  in  standard  tables  together  with  the  spectral  Irrad- 
iance  of  the  three  C.I.  E.  standard  sources  (A,B,C). 

The  standard  chromaticity  diagram  (Fig.  1)  provides  a two-dimensional 
characterization  of  a color  stimulus.  The  "chromaticity"  coordinates 
are  given  by: 

X = X ; y = Y 

X Y Z X Y Z 

On  the  diagram  the  colors  obtained  by  mixing  two  others  will  lie  on  a 
straight  line  connecting  the  two  points.  The  curved  boundary  line  is 
the  locus  of  the  pure  spectral  colors,  and  the  area  within  this  boundary 
contains  all  the  colors  that  can  be  made  with  lights,  with  white  at  the 
center. 

The  simplest  theory  of  the  mechanism  of  color  vision,  proposed  by 
Young  and  Helmholtz,  postulates  that  there  are  exactly  three  different 
receptor  pigments  in  the  eye,  which  respond  maximally  to,  say,  red 
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green,  and  blue  wavelengths.  These  three  channels  of  information 
are  then  processed  somehow  by  the  visual  system  so  that  colors  can  be 
perceived. 

The  laws  of  color  mixing  can  be  demonstrated  to  be  a consequence  of 
this  proposition,  but  since  we  can  make  any  kind  of  linear  transforma- 
tion of  the  color  coordinates  that  we  want,  we  cannot  determine  uniquely 
what  the  three  receptor  absorption  curves  are  simply  on  the  basis  of 
color  mixing  experiments. 

COLOR  RECEPTORS 

For  several  decades  a number  of  workers  have  focused  their  attention 
on  the  problem  of  determining  the  nature  of  the  visual  photopigments. 

The  most  definitive  data  has  come  from  the  microspectrophotometry 
of  Marks,  Dobelle,  and  MacNichol,  in  1964,  and  similar  work  by 
Brown  and  Wald,  1964.  Marks  et.  al.  (1964)  present  data  from  measure- 
ments taken  on  ten  primate  cones  and  two  human  cones  from  the  para- 
foveal region.  Their  procedure  included  the  use  of  an  "end-on"  pre- 
paration of  retinal  receptors,  so  that  the  measuring  light  beam  passed 
through  the  receptor  in  the  normal  direction.  They  simultaneously 
passed  a reference  beam  through  a "blank"  spot  in  the  preparation, 
and  the  ratio  of  the  intensity  of  the  two  beams  was  measured  auto- 
matically over  the  range  of  390  nm  to  67  0 nm. 

An  overview  of  the  results,  including  primate  cones,  reveals  three 
major  classes  of  receptors  with  maximum  absorption  at  about  445, 

535,  and  57  0 nm.  One  of  the  human  cones  showed  maximum  absorp- 
tion at  457  nm,  indicating  that  it  was  a mediator  for  blue  sensation. 

The  other  human  cone  tested  by  Marks,  et.  al.  (1964)  was  probably 
a red  receptor,  with  its  peak  absorption  at  57  5 nm. 
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The  absorption  difference  spectra,  from  the  majority  of  the  "red" 
cones  showed  a secondary  hump  or  peak  near  555  nm  suggesting  the 
possibility  that  single  "red"  receptors  may  contain  red  and  green 
pigments,  coexisting  in  a single  cone. 

Brown  and  Wald  (1964)  measured  the  difference  spectra  of  4 human 
cones  using  a similar  "end-on"  preparation  from  the  parafoveal 
region  of  the  retina.  Their  technique  differed  from  Marks,  et.al.  (1964) 
in  that  they  took  a difference  spectrum,  i.e.  , they  obtained  the  absorp- 
tion spectrum,  then  bleached  the  receptors  with  a bright  flash  of 
light  and  obtained  a bleached  absorption  spectrum.  The  difference 
between  the  unbleached  and  bleached  absorption  spectra  is  the  difference 
spectra.  Brown  and  Wald's  data  represent  one  blue  receptor,  two  green- 
receptors,  and  one  red  receptor.  The  peaks  in  these  curves  fell  at 
about  450,  525,  and  555  nm.  The  blue-receptor  peak  of  450  nm  agrees 
well  with  Marks  et.al.  blue  receptor  which  peaked  at  457  nm,  but  the 
pree  and  red-receptor  peaks  seem  to  occur  at  somewhat  shorter  wave- 
lenghts  than  was  found  earlier  by  Brown  and  Wald  when  relatively 
large  patches  of  human  fovea  were  measured. 

TETRACHROMATICITY 

It  is  often  stated  that  there  are  four  "psychologically  unique"  colors; 
red,  yellow,  green  and  blue.  These  lead  to  the  so-called  opponent- 
color  theories,  such  as  those  proposed  by  Hurvich  and  Jameson  (1955). 
The  blue,  green  and  yellow  primaries  are  also  special  in  that  they 
do  not  undergo  a Bezold-Brucke  shift;  that  is,  their  hue  is  invariant 
with  changes  in  luminance.  Psychometric  experiments  show  red  and 
green  as  antagonistic  colors;  their  mixture,  yellow,  is  not  perceived 
by  subjects  as  either  red  nor  green.  Conversely,  blue  and  yellow 
lights  mix  to  form  white.  Hurvich  and  Janeson  (1955)  have  worked  out 
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a formal  mathematical  theory,  based  on  the  opponent -color  hypothesis 
with  which  they  are  able  to  match  predicted  and  experimental  data 
quite  well  for  various  phenomena,  including  the  Bezold-Brucke  shift, 
wavelength  discrimination,  spectral  saturation,  and  color  mixture. 
These  results  are  not  surprising  since  their  chromatic  response 
functions  are  linear  transformations  of  the  CIE  mixture  curves  with 
the  addition  of  four  "variable  constants"  which  are  manipulated  for 
best  fit. 

Color  theory  is  still  evolving,  and  it  was  felt  that  this  was  an  area  in 
which  evoked  response  studies  could  make  an  important  contribution. 
This  choice  provided  an  extremely  fortunate  collaboration  with 
Dr.  C.T.  White  at  NEL.  Indeed,  his  preliminary  data  indicated  that 
a paradigm  of  this  type  would  be  a very  rich  vein  for  our  purpose. 

A sample  of  that  preliminary  datais  shown  in  Fig.  2.  There, 
average  evoked  responses  are  shown  to  red,  green  and  blue  flashes 
presented  against  yellow  backgrounds.  Background  level  on  the  one 
hand  and  stimulus  intensities  in  each  color  on  the  other  hand,  were 
varied  separately.  From  the  ensemble  of  these  experiments  obtained 
under  different  background  colors  and  intensities,  tentative  candidate 
codes  emerged  that  accounted  for  the  three  hypothetical  processes 
(at  retinal  level)  corresponding  to  the  red-green-blue  trichromatic 
theory  (Fig.  3).  The  yellow  background  was  thought  to  produce 
particularly  representative  evoked  responses  because  of  the  red- 
green  depression  caused  by  retinal  adaption  to  yellow  light.  This 
allows  the  blue  component,  which  had  been  elusive  before,  to  appear 
much  more  clearly.  These  were  startling  results  since,  while 
there  is  good  evidence  for  a trichromatic  absorbing  structure  in 
the  fovea  (Fig.  4)  there  was  no  indication  that  these  processes 
would  be  reflected  into  scalp  potentials  (collected  in  the  occipital 
area).  This  finding  also  brings  as  a next  question,  the  potential 
role  of  the  opponent  color  theory  developed  from  psychophysical 
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...  Curves  showing  dependence  of  magnitudes 
of  the  three  basic  retinal  processes  upon  wavelengths 
of  spectral  lights,  obtained  from  fovea  of  normal  tri- 
chromat. 'MOTOKAWA,  1919. 


experiments. 

Topological  Dimensionality  studies  of  this  early  data  showed  that 
red  and  green  components,  although  very  close  in  timing,  did  behave 
in  a clearly  distinct  and  qua  si -orthogonal  manner  with  respect  to 
each  other  when  background  or  stimulus  intensities  were  changed 
(Fig.  5).  The  response  of  the  dark  adapted  eye,  a limit  case  in 
the  background  level  series,  clearly  detached  itself  from  the  cluster 
in  all  cases.  A fraction  of  this  preliminary  data  was  obtained  at  NEL 
without  the  benefit  of  automatic  data  acquisition  and  digitizing 
systems.  Thus,  single  epochs  had  been  lost  and  only  graphic  aver- 
ages were  available.  Further  analysis  of  these  averages  has  been  pursued  with 
borrowed  graph  pen  to  digitize  plotted  data.  Since  then,  the  experi- 
ments have  been  conducted  at  the  Brain  Computer  Interface  Labora- 
tory and  have  been  or  are  being  analyzed  through  the  BCI  software 
package  without  resorting  to  averaging  (except  for  the  final  displays). 
Twelve  subjects  have  been  evaluated  on  the  basis  of  evoked  response 
quality  and  relative  absence  of  muscle  artifacts  and  fivehayo  been  retained 
retained.  The  BCI  experiments  are  focused  on  single  epochs  (aver- 
ages are  actually  irrelevant  to  the  Biocybernetic  goals.) 

Meanwhile,  on-line  experiments  have  been  done ' to  further  probe 
the  fine  structure  of  the  EEG  epochs  and  search  for  better  primitives 
than  amplitude  samples.  Particularly  intriguing  was  the  possibility 
that  by  adequate  representation  of  the  data,  these  primitives  could  be 
made  explicit.  Various  representations  were  attempted,  among?  which 
the  replacement  of  the  single  evoked  response  by  the  point  processes 
formed  by  the  instants  at  which  a single  epoch  was  displaying  a peak 
in  the  positive  or  negative  direction  (Fig.  6).  In  other  words,  the 
series  of  zero  crossing  instants  for  the  first  derivative  of  the  ampli- 
tude function  together  with  the  sign  of  the  second  derivative.  Cumulative 
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histograms  of  those  micro-events  were  obtained  and  compared  to 
the  averages  (Fig.  7).  These  histograms  were  found  to  be  as 
stimulus -specific  as  are  the  averages.  While  the  relationship 
between  the  histograms  and  the  underlying  slow  wave,  can  be 
explained  to  some  extent  by  the  procedure  itself  (i.e.  the  ’ bunching" 
of  the  peaks  by  the  slow  wave),  it  has  already  been  established  that  the 
stability  and  sharpness  of  some  of  the  peaks  could  not  be  explained 
if  the  small  peaks  were  just  high  frequency  noise.  Rather  it  appears 
that  the  slow  wave  conjugates  with  short  events  which  may  corre- 
spond to  the  firing  of  specific  nuclei.  This  raises  great  hope  that 
a combination  of  short  event  detection  with  amplitudes  at  selected 
times  will  provide  encoding  of  the  evoked  response  in  a particularly 
economical  manner  in  man-machine  communication.  Already  at 
this  time,  even  without  the  benefit  of  amplitude  information  the  peak 
"message"  has  equivalent  clarity  or  discriminability , when  compared 
to  the  peak  histogram,  as  the  single  epoch  has  with  respect  to  the 
average. 

A next  move  is  underway  to  investigate  the  possibility  that  the  peak 
set  or  some  variation  of  it  contains  identifiable  moment-to-moment 
information.  Operant  conditioning  techniques  are  expected  to  provide 
some  of  the  answers  and  various  operant  paradigms  are  being  evaluated. 
Present  results  are  still  at  a level  of  casual  observations.  Fig. 
for  example , shows  a run  in  which  the  reward  (a  conditioning)  was 
attached  to  the  presence  of  small  positive  peaks  in  the  vicinity  of 
300  milliseconds.  The  growth  of  the  amplitude  of  this  component  in 
the  average  is  clearly  seen  in  the  graph. 
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DATA  ANALYSIS 


Phase  I 


Data  from  one  subject  whose  Evoked  Potentials  (EVPOTS)  were  extensively 
studied  is  presented  next.  This  experiment  consisted  of  48  presentations  of 
a flash  of  colored  light  of  either  red,  green  or  blue  hue,  over  a yellow  back- 
ground. The  brain  electrical  activity  was  sampled  for  5 seconds  prior  to  and 
320  milliseconds  following  each  flash,  from  four  electrode  locations.  A 
description  of  the  experiment  follows: 


Exp  ID: 
Stimuli: 


No.  of  Epochs 
per  group: 
Background: 
Sampling  Rate: 
Channels  Recorded: 


Three  Color  Experiment 
Flashes  of  red,  green  and  blue.  The  flashes 
were  non-random,  i.e.  three  sequences  of 
color  flashes  were  used  to  record  three 
batches  of  evoked  responses 


48 

Yellow 
320  msec. 
Channel  1 
Channel  2 
Channel  3 
Channel  4 


Frontal/Ears 

Vertex/Ears 

Occipital/Ears 


,+A 


PZ 

o,/a;+a 


1 h2 


Cz/Aj+A, 


(Eye)  Inferior  Orbital/E?ars 


This  experiment  is  not  intended  to  represent  the  definitive  color  evoked 
lesponse  study;  rather  it  is  a vehicle  for  the  evaluation  of  evoked  response 
analysis  techniques.  The  data  have  been  analysed  using  stepwise  discrim- 
inant analysis  (SWDA)  with  and  without  pre-processing  by  the  Orthogonal 
Transformation. 


Some  representative  evoked  response  plots  are  shown  in  the  following  figures. 
Fig.  9 in  the  plot  of  Red,  Green  and  Blue  EVPOTS  from  Channel  1,  Fig.  10 
is  from  Channel  2,  and  Fig.  11  is  from  Channel  3.  These  plots  are  all  aver- 
ages obtained  from  48  epochs  of  each  of  the  three  colors.  The  ordinate 
is  in  microvolts,  positive  up,  the  abscissa  is  time;  there  are  80  samples 
spaced  4 msec,  apart.  The  flash  occurred  at  time  point  1. 


{Frontal 1 Ears) 


Clearly,  the  average  EVPOTS  evoked  by  the  different  colors  exhibit  differences 
but  we  need  a measure  of  how  different  they  are,  with  an  appreciation  of 
the  variability  from  trial  to  trial.  The  stepwise  discriminant  analysis  (SWDA) 
program  does  provide  an  objective  measure  of  the  differences  between 
these  multivariate  distributions,  based  on  the  variability  and  the  regularity 
in  the  data. 

The  program  computes  a set  of  linear  classification  functions  by  choosing  a 
subset  of  the  original  80  variables  in  a stepwise  manner.  The  variable  entered 
at  each  step  is  selected  according  to  one  of  four  possible  criteria;  typically 
that  variable  will  be  chosen  which  has  the  largest  F-to-enter.  The  F statistic 
which  is  computed  at  each  step  gives  us  an  indication  of  the  differences  between 
groups  relative  to  the  differences  within  groups. 

Typically,  the  SWDA  program  is  told  to  select  the  10  best  variables,  then  it 
computes  the  coefficients  of  these  variables  which  produces  the  best  possible 
discrimination  between  groups,  and  prints  a classification  matrix.  Each  case 
(epoch)  is  evaluated  using  the  discriminant  function  (D.  F.)  developed  and  its 
posterior  probability  of  belonging  to  each  group  (color)  is  computed  as  well  as 
the  square  of  the  Mahalanobis  distance  from  each  group.  The  program  also 
computes  the  coefficients  for  canonical  variables  and  plots  the  first  two  can- 
onical variables  to  give  an  optimal  two-dimensional  picture  of  the  separation 
of  groups. 

Other  measures  generated  by  SWDA  include  group  means  and  standard  deviations, 
within  groups  covariance  matrix,  the  within  groups  correlation  matrix;  and  for 
each  step,  the  variables  included  and  F-to-remove,  variables  not  included  and 
F-to-enter,  Wilk's  Lambda  (or  U statistic)  and  approximate  F statistic  to  test 
equality  of  gtoup  means,  and  the  matrix  of  F statistics  to  test  the  equality  of 
means  between  each  pair  of  groups. 
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A subset  of  this  wealth  of  information  about  the  data  has  been  tabulated 
and  will  be  referred  to  in  the  following  paragraphs,  along  with  comments 
about  it  significance. 

Referring  to  the  map  distances  in  Table  1,  it  is  seen  that  blue  and  red  are 
the  closest  pair  in  Ch.  3 (over  visual  cortex),  while  blue  and  green  are 
closest  for  Ch.  4. 

This  finding,  if  further  validated,  may  lead  to  a demonstration  of  an 
objective  measure  of  cortical  processing  of  hue  information  in  intact  man. 

This  finding  could  be  predicted  by  pointing  out  that  blue  and  red  which  are 
at  the  two  extremes  of  the  spectrum  are  also  least  similar  at  the  eye  level  (Ch.  4) 
which  in  part  reflects  the  EOG  (Electrooculogram).  By  contrast,  the  nearness 
of  blue  and  red  over  visual  cortex  (Ch.  3)  matches  the  data  in  the  CIE  chro- 
maticity  diagram  (Fig.  1)  which  was  constructed  from  psychophysical  experi- 
ments, thus  probably  reflecting  processes  occurring  in  the  visual  cortex. 

Also  from  Table  1,  referring  to  the  10  step  classification  matrices,  it 
can  be  seen  that  SWDA  did  a very  good  job  of  classifying  each  epoch  into  its 
proper  group  using  a D.  F.  based  on  only  10  sample  variables.  Ch.  3 gave 
the  most  accurate  classification,  with  an  error  rate  of  7.6  when  operating 
with  three  groups  with  a misclassification  rate  of  8.  9 . The  classification 
matrix  also  provides  distance  information  in  the  form  of  confusibility. 

Looking  at  the  10  step  matrix  for  3 colors,  for  Ch.  3 (Table  1),  it  is  seen 
that  green  epochs  are  always  classified  as  belonging  to  the  green  group,  48 
out  of  48,  while  44  red  epochs  were  correctly  classified  as  red,  none  as  green, 
and  4 were  classified  as  blue.  Also,  7 epochs  of  blue  were  incorrectly  class- 
ified as  red,  none  as  green,  and  40  were  correctly  placed  in  the  blue  column. 
These  results  imply  that  red  and  blue  are  closer  than  either  red  and  green  or 
blue  and  green  in  these  data,  since  they  are  more  often  confused,  i.e.  incorrectly 
classified  as  each  other  by  the  program. 
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A.  Channel-Wise  Comparisons  (Table  1) 

1)  As  expected,  with  the  data,  Channel  3 gave  the  best  discrimina- 
tion followed  by  Channels  4,  2 and  1 respectively. 

2)  Channels  1 and  2 chose  similar  initial  variables  to  discriminate 
between  the  different  color  groups. 

0 

3)  The  time  points  around  26,  29,  30,  33,  41,  45,  61,  and  69 
seemed  to  predominate  the  classification  functions  computed 
for  the  various  channels  from  the  raw  data. 

0 

B.  Color-Wise  Comparisons 

l.  Red  (Table  3 run  A-9,  Table  4 runs  A-8,  A-12a 

1)  Red  color  data  is  easily  separable. 

0 

2)  Using  the  Channel  3 runs  (10,  16,  18)  we  see  that  time  points 
around:  21,  29,  33,  36,  46  predominate. 

3)  We  also  note  that  time  point _29  occurs  in  all  the  three  runs 

indicating  that  time  point  29  definitely  is  associated  with  a J) 

red  color  component. 

4)  Overall,  time  points  21,  29  and  33  can  be  attributed  to  the 
color  red. 

O 

n.  Green  (Table  3 run  A-12,  Table  4 runs  A-12a,  A-7) 

1)  Green  color  seems  to  be  easily  distinguishable  by  SWDA. 

2)  Using  Channel  3 runs  (12,  18,  20)  we  see  that  time  points 
around:  10,  20,  28,  41,  45,  64,  73  seem  to  dominate. 

3)  Overall,  time  points  28,  45  could  be  attributed  to  the  color  green. 

m.  Blue  (Table  3 run  A-ll,  Table  4 runs  A-8,  A-7) 

1)  Using  the  Channel  3 runs  we  find  that  time  points  around; 

18,  21,  22,  24,  33,  37,  41,  79  seen  to  dominate. 

2)  Time  point  22,  38  gives  us  a hint  of  the  color  blue. 

O 
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C.  Inter  Color  Qualifications 


I.  Red  vs  Green  (Tables  3 and  4 runs  A-9,  A-12,  A-12a) 

Using  Channel  3 we  see  that  time  points  around;  10,  20,  29, 

33,  34,  36,  40,  64,  73  dominate. 

n.  Green  vs  Blue  (Tables  3 and  4 runs  A-12,  A-ll,  A-7) 

Using  Channel  3 we  see  that  time  points  around:  10,  20,  33, 

34,  37,  41,  45,  64  dominate. 

in.  Red  vs  Blue  (Tables  3 and  4 runs  A-9,  A-ll,  A-8) 

Using  Channel  3 we  find  that  time  points  around:  20,  26,  29 
33,  34  dominate. 

To  answer  the  Question  as  to  how  good  the  classification  functions  produced 
by  SWDA  are,  a special  run  was  made.  The  data  used  were  from  Chennel  3 
(raw).  The  Red,  Green  and  Blue  groups  were  each  split  in  two  groups,  and 
the  classification  functions  produced  by  the  program  based  on  the  first  half 
of  the  data  were  used  to  classify  the  second  half.  The  classification  matrices 
at  steps  1,  5 and  10  are  given  in  Table  5. 

We  see  that  the  classification  functions  produced  by  the  program  provide  a 
qood  discrimination  (86  vs  93  overall)  on  color  evoked  responses.  This 
provided  encouragement  toward  the  use  of  the  classification  functions 

on-line  in  subsequent  experiments. 
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COMPARISONS  TABLE  #2 


COMPARISONS  TABLE  #3 
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A pictorial  display  of  the  classification  and  separation  ability  of  SWDA 
is  provided  in  the  form  of  the  now  classic  canonical  separation  map. 

Each  epoch  is  located  in  the  two-dimensional  space  and  labeled  with 
the  first  letter  of  its  group  name.  An  asterisk  is  placed  at  the  center 
of  "mass"  of  each  group,  and  overlaps  are  indicated  by  a $. 

Figures  21  thru  24  are  representative  canonical  maps  which  allow 
further  comparisons  of  the  effect  of  channel,  color,  orthogonalization 
and  number  of  steps  on  the  separation. 

CONCLUSIONS  : Phase  I 

The  use  of  stepwise  discriminant  analysis  is  valuable  for  the  classifica- 
tion of  evoked  potential  data.  The  program  efficiently  selects  that  sub- 
set of  samples  which  can  accurately  characterize  the  data  in  an  optimum 
manner.  The  samples  which  are  selected  can  be  associated  with  the 
"components"  of  the  evoked  potential  waveform  to  which  most  researchers 
have  paid  attention  in  previous  work.  Furthermore,  the  discriminant 
function  developed  by  SWDA  can  be  used  for  efficient  and  highly  accurate 
on-line  classification  of  single  epochs  of  evoked  potential  data,  for  use 
in  operant  conditioning  experiments. 

The  limitations  of  SWDA  arise  from  the  nature  of  the  model  which  is 
presumed  to  underly  the  generation  of  evoked  potentials,  that  is  one 
in  which  a deterministic  signal  is  added  to  random  noise.  Other  models 
of  the  EVPOT  can  be  entertained,  for  which  other  methods  of  analysis 
would  be  appropriate.  The  detection  of  "peaks"  (points  at  which 
dv/dt  = 0)  is  one  such  alternative  approach  which  we  are  actively 
pursuing.  At  the  present  state  of  the  art,  it  would  seem  that  a combin- 
ation of  peak  detection  techniques  and  amplitude  measures  could  provide 
the  most  useful  and  elegant  way  of  representing  the  EVPOT  "signature". 
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data  analysis 


PHASE  2 

The  search  for  an  appropriate  methodology  was  centered  on  what  appeared  to  be 
"meaningful"  questions  relative  to  color  evoked  responses.  A summary  of  those 
questions  is  given  in  Table  8 along  with  the  run  codes  of  the  computer  runs 
which  are  summarized  later  In  the  report. 

The  sheets  in  the  Appendix  are  the  condensations  of  the  results  of  Stepwise 
Discriminant  Analysis  runs  on  data  provided  by  basically  two  experiment  groups. 
The  two  experiments  are  described  below: 


Experiment  1:  Ten  Color  Experiment 

Subject:  CP 

Background:  250  Ft  Lamberts-Yel low 

EEG  Channels:  1-70  Hz 

Linked  Ear  Reference: 

/ 

Top  View 


A 3 A ^2 


Sensitivity:  2uv/nm 


Sampling  Rate:  4 msec. 

Total  Epoch  Duration:  320  msec. 

Stimuli  were  non-random,  i.e.  64  epochs  for  each  color  stimulus 
were  taken  i ;i  sequence.  For  order  of  stimuli  see  Table  6 

No.  of  Epochs  per  color:  64 

No.  of  Sample  Points/epoch:  80 


No.  of  Channels:  5 


Channel  1:  Top  Occipital 

Channel  2:  Central  Occipi-a! 

Channel  3:  Bottom  Occipital 

Channel  4:  Left  Occipital 

Channel  5:  Right  Occipital 


3 cm  ant  to  0 /A 
07/A 

3 cm  post,  to  0 /A 
0,/A 
O^/A 
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In  phase  2,  only  the  data  belonging  to  the  briqhter  stimuli,  i.e.  at  0.0 
loq  unit  attenuation  (namely  data  sets  1 through  A)  were  processed. 

This  experiment  was  performed  mainly  to  see  the  effects  of  spanning  the 
whole  spectrum  thus  obtaining  a base  for  analysis.  The  approach  was 
also  a search  for  components  associated  with  the  three  basic  processes, 
namely,  red,  green  and  blue  and  for  investigating  of  the  wrap-around  effects 
of  the  spectrum  if  present  in  the  E.R.  The  experiment  gave  evoked  responses 
on  a larqe  variety  of  colors  thus  providing  numerous  combinations  to  search 
for  component  codes. 

Experiment  2:  Opponent  Color  Experiment 
Subject:  CP 

Sensitivity:  2uv/mm  (Except  for  Purple  which  was  3 uv/mm) 

Sampling  Rate:  4 msec. 

Total  Epoch  Duration:  320  msec 

Stimuli  were  non-random,  i.e.  50  epochs  1 

for  each  color  stimulus  were  taken  In  sequence.  For  order  of 
stimuli  see  Table  7 TX. 


No.  of  Epochs  per  color:  50 
No.  of  Sample  Points/Epoch:  80 
No.  of  Channels:  5 

Channel  1 : 
Channel  2: 
Channel  3: 
Channel  4: 
Channel  5: 


Top  Occipital  3 cm  ant.  to  0 /A 

Central  Occipital  , 0 /A 

Right  Occipital  02/A 

Left  Occipital  O./A 

Eye 
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TABLE  7 


OPPONENT  COLOR  EXPERIMENT 


Exp.  ID 

Color 

Attenuation 

502001 

Blue 

2.5 

502002 

(465 

nm) 

2.0 

502003 

1.5 

502004 

1.0 

502005 

0.5 

502006 

0.3 

502007 

> 

/ 

0.0 

502008 

Yellow 

2.5 

502009 

(575nm) 

2.0 

502010 

1.5 

50201 1 

1.0 

502012 

0.5 

502013 

0.3 

502014 

\ 

r 

0.0 

502015 

Green 

2.5 

502016 

(515nm) 

2.0 

502017 

1.5 

502018 

1.0 

502019 

0.5 

502020 

0.3 

502021 

> 

* 

0.0 

502022 

P 

ed 

2.5 

502023 

(62C 

nm) 

2.0 

502024 

1.5 

502025 

1.0 

502026 

0.5 

502027 

0.2 

502028 

% 

r 

0.0 

510001 

Purple 

0.0  to  6 

510002 

0.3 

510003 

0.5 

510004 

1.0 

510005 

1.5 

510006 

> 

¥ 

2.0 

510007 

Blue  Green 

2.5 

510008 

(490nm) 

2.0 

510009 

1.5 

510010 

1.0 

510011 

0.5 

510012 

0.3 

510013 

> 

f 

0.0 

510014 

White 

1.8  (1+.5+.3) 

510015 

1.5  (1+.5+.1) 

510016 

1.5  (1+.5) 

510017 

1.9  (1+.5+.3+.1) 

510018 

> 

V 

1.7  (1+.3+.3+.1 ) 

Sensitivity 


2uv/mm 


* 

3uv/mm 


2uv/nui 


♦ 
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TABLE  8 EXPLORATORY  DATA  ANALYSIS  USING  DATA  OBTAINED  FROM  THE 
TEN  AND  OPPONENT  COLOR  EXPERIMENTS 


QUESTIONS  POSED 


1.  Ouantlf Icatlon  of  dl scrim-inability  of  different 
datum  present  In  the  early  period  of  the  epoch. 

2.  Usaqe  of  multiple  channels  In  the  training  and 
testlnq  of  DF's. 

3.  Usaqe  of  a low  band  pass  hlqh  qaln  filter  prior 
to  training  the  DF. 

4.  DF  trained  on  data  from  one  session  and  tested 
on  data  from  another. 

5.  DF  trained  on  data  of  one  Intensity  and  applied 
on  data  of  hlqher  or  lower  Intensity. 

6.  DF  trained  on  averages  and  applied  on  raw  data. 

7.  DF  trained  on  single  ER  and  applied  on  slnqle 
ER. 

8.  To  study  the  variation  in  ER  as  the  whole 
spectrum  Is  spanned. 

9.  To  study  the  closeness  of  Greefj  Green-Yellow 
and  Yellow-Green  In  comparison  to  the  extremes 
of  the  spectrum. 

10.  To  Investigate  wrap-around  effects  of  the 
spectrum. 


11.  To  study  the  performance  of  DF  trained  on  Red, 
Green  and  Blue  ER's  and  applied  on  ER's 
recorded  from  flashes  of  others  colors. 

12.  Investigation  of  the  dlscrimlnabilitv  of  ER's 
obtained  with  opponent  color  pairs 

13.  To  study  the  dlscrlmlnabillty  of  opponent  colors 
of  Red,  Green  and  Blue,  i.e.  Blue-Green,  Purple 
and  Yellow  respectively. 

14.  In  search  for  component  codes  by  discriminating 
aqalnst  a constant  DC  response. 

15.  In  search  for  component  codes  by  discriminating 
against  ER's  recorded  for  White  color 


RUNS  MADE  (see  Appendix) 

B-l 

B-l,  C-l,  C-2 
C-4,  C-5 
D-4 

D-13,  D-14,  D-15 

C-7 , C-9,  C-10,  C-l 1 , C-l la,  C-8 

C-l,  C-2,  C-3,  C-4,  C-8,  C-l 2, 
D-4,  D-6,  D-7,  D-8,  D-9,  D-10, 
D-ll,  D-12,  D-19,  D-21 , D-22,  D 

B-3 , B-4,  B-5,  B-6,  B-7 
C-3,  C-4 


B-3,  B-4,  B-5,  B-6,  B-7,  B-8, 
B-9,  C-l,  C-2,  C-3,  C-4,  C-5 


B-3,  B-4,  B-5,  B 5,  B-7,  B-8, 
B-9 


5 

a 

( 

| 


-23 


D-ll,  D-12,  D-22 


D-19,  D-21 

C-l 3,  C-l 3a , C-l 4 , C-l 6,  C-18 
C-l 9,  C-20 

D-6,  D-7,  D-8,  D-9,  D-10 


'mb 


QUESTIONS  POSED 


RUNS  MADE  (see  Appendix) 


16. 

To  investigate  the  possibility  of  expressing 
ER's  due  to  one  color  In  terms  of  two  others 
in  the  case  of  tripole  of  Red,  Green  and  Blue. 

c-in,  C-ll,  C-lla 

17. 

To  evaluate  the  window  estimated  by  runs 
associated  with  Red,  Green  and  Blue  colors. 

D- 23 

18. 

To  investigate  the  effect  on  the  predictive 

power  of  DF  when  identified  outliers  are 
excluded  from  the  computation  of  DF 

D-l , D-2,  D-3 
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The  experiment  was  run  using  seven  different  intensity  levels  (from  0 log 
units  attenuation  to  2.5)  for  each  color.  A few  runs  were  made  to  investi- 
gate the  "goodness"  of  a discriminant  function  developed  on  a training  set 
at  one  intensity  and  applying  It  to  the  testing  set  of  another.  Besides 
these,  all  other  exploration  was  restricted  to  the  data  sets  involving  the 
brightest  intensity  in  order  to  have  proper  comparisons  with  the  data 
obtained  from  the  ten  color  experiment. 

We  shall  now  go  into  the  details  of  the  runs  made  and  questions  asked  as 
listed  in  Table  8 . The  runs  referred  to  are  summarized  later. 

QUESTION  1;  R jn  B-l 

To  investigate  the  amount  of  information  contained  in  the  early  period 
j of  the  epoch  (i.e.  right  after  the  stimulus)  reflecting  the  separability 

of  the  different  groups. 

A stepwise  discriminant  analysis  was  performed  on  three  colors.  The  discrim 
inant  functions  (DF)  were  obtained  on  the  first  10  sample  points  only,  i.e. 
the  first  40  msec,  of  the  epoch. 

As  we  see  from  the  run,  the  discriminability  of  the  data  was  estimated  at 
47.92%  and  the  predictive  power  of  the  OF  was  measured  to  be  only  34.38%. 

} 

The  poor  performance  is  at  chance  level.  Clearly  there  is  very  little  if 
any  information  relevant  to  the  color  groups  contained  in  the  early  period 
of  the  evoked  response. 

if 

This  confirmed  data  from  other  sources  and  on  this  basis,  most  of  the  subse- 
quent runs  were  made  ignoring  the  first  10  sample  points. 
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QUESTION  2:  Runs  B-l.  C-1.  C-2 

To  Investigate  the  possibility  of  determining  a better  performing  DF 
by  using  two  central  occipital  channels  (rather  then  one)  to  train  the 
OF.  (The  electrodes  are  one  Inch  apart.) 

The  first  run,  B-l,  was  made  using  only  the  early  portion  of  the  epoch 
and  aid  not  give  any  useful  Information.  However,  C-1  was  run  using 
the  Top  and  Central  Occipital /Ears  channels  whereas  C-2  was  run  using 
only  the  Top  Occipital/Ears  channel.  In  C-1,  all  the  odd  sample  points 
were  from  CHI  (Top  Occipital /Ears)  and  the  even  ones  were  from  CH2 
(Central  Occipital /Ears). 

As  seen  in  Table  9 below  there  was  no  advantage  demonstrated  by  using 
the  two  channels  In  training  and  testing  the  DF.  Rather  the  performance 
remained  unchanged.  This  may  be  due  to  the  close  proximity  of  the 
electrodes  In  which  case  It  would  suggest  a lower  bound  for  Interelectrode 
distance.  Again,  this  is  In  agreement  with  general  working  assumptions 
In  the  field. 

TABLE  9 USAGE  OF  REDUNDANT  CHANNELS  IN  THE  TRAINING  & TESTING  OF  DF'S 


C 


O 


O 


Qualified 

CHI  + CH2 

CHI  Only 

1. 

Discrlmlnabllity  of  data 

93.75% 

93.75% 

2. 

Predictive  Power  of  DF 

84.88% 

84.88% 

3. 

Initial  Approximate  F 

31.38367 

31.67860 

4. 

Variables  Selected 

CHI: 

63,49,21,  23 

CHI:  63,  31,  52, 

45 

CH2: 

38,30,28,  44, 

48,  24,  21  , 

28, 

46. 

66,  55. 

O 
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QUESTION  3;  Runs  C-4,  C-5 

To  study  the  effects  of  filtering  on  the  discriminate ity  of  data  and 
subsequent  performance  of  the  DF. 

A double  pole,  low  pass  digital  filter  with  the  ability  to  zero  mean  and 
time  varying  "time  constants"  from  Glassman's  Orthogonal  Derivation  Program 
was  programmed  to  pre-process  the  data  before  the  Stepwise  Discriminant 
Analysis. 

The  data  consisted  of  a training  set  of  32  epochs  each  of  Violet,  Green, 
Green-Yellow,  Yellow-Green  and  Red  as  well  as  a testing  set  of  similar 
contents. 

Two  runs  were  made  of  the  Stepwise  Discriminant  Analysis.  In  one  case, 
the  data  had  been  pre-processed  by  the  filter  and  in  the  other  the  raw 
data  was  analyzed. 

As  a result,  both  the  discriminability  of  the  data  (from  66%  to  70%)  and 
the  predictive  power  of  DF  (from  58%  to  64%)  improved  suggesting  that 
filtering  of  the  data  may  be  beneficial.  Following  are  the  different 
variables  (time  points  in  msec/4)  in  ascending  order  for  the  two  runs: 

Non-Filtered:  24,  28,  30,  33,  37,  42,  46,  49,  54,  56; 

Filtered:  22,  24,  25,  28,  33,  38,  45,  49,  55,  66; 

Time-wise  the  variables  selected  fall  in  the  same  general  neighborhoods. 

Further  work  on  applicability  of  the  filter  and  optimization  of  the 
parameters  was  planned.  The  lowpass  value  for  the  filter  was  adjusted 
to  match  that  of  the  EEG  polygraph.  The  improvement  is  therefore 
attributed  to  a sharper  drop-off  of  the  transfer  function  at  high  fre- 
quencies. 
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QUESTION  4;  Run  D-4 


To  study  effects  of  training  the  DF  on  data  of  one  session  and  applying 
it  to  that  of  another  session  with  the  same  subject. 


In  this  run  a DF  was  developed  on  32  epochs  each  of  Red,  Green,  Blue 
and  Yellow  from  the  Ten  Color  Experiment  and  applied  on  50  epochs  each 
of  Red,  Green,  Blue  and  Yellow  of  the  Opponent  Color  Experiment.  This 
was  done  to  investigate  the  validity  of  DF  between  sessions.  At  this 
the  results  were  not  very  encouraging.  The  discriminability  of  data 
in  the  Ten  Color  Experiment  was  only  70.31%  and  the  predictive  power  of  the 
DF  system  obtained  when  applied  on  the  Opponent  Color  data  fell  to  40.50% 


Apparently  the  DF  obtained  was  negatively  affected  by  the  introduction 
of  Yellow  data,  as  without  it,  data  from  the  same  experiment  using  Red, 
Green  and  Blue  had  yielded  discriminability  of  91.67%  (Run  D-l). 
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QUESTION  5;  Runs  D-13,  D-14,  D-15 
To  study  effect  of  intensity  on  DF. 

Three  runs  were  made  using  different  combinations  of  intensity: 

D-13  DF  trained  on  high  intensity  (0.0) 

and  applied  on  low  intensity  (0.3,  0.5) 

D-14  DF  trained  on  low  intensity  (0.3) 

and  applied  on  higher  intensity  (0.0)  and 
on  lower  intensity  (0.5) 

D-15  DF  trained  on  low  intensity  (0.5) 

and  applied  on  high  intensity  (0.3)  and 
higher  intensity  (0.0) 


TABLE  10  QUANTIFICATION  (percentages)  OF  DISCRI MINABIL ITY  OF  DATA 
(diagonal)  AND  PREDICTIVE  POWER  OF  DF  (non-diagonal) 


Training/ testing 
Intensity  / Intensity 

0.0 

0.3 

0.5 

(Brightest)O.O 

86 

55 

30 

0.3 

46 

83 

49 

0.5 

31 

51 

57 

Therefore  1)  The  discriminability  of  data  degrades  as  the  intensity  of 
stimulus  decreases  as  would  be  predicted  by  considerations  of  signal-to- 
noise  ratio.  2)  The  discriminant  function  are  intensity-bound  and  do 
not  generalize.  This  is  not  surprising  since  the  DF  primitives  are 
time-samples.  Intensity  changes  are  likely  to  bring  about  significant 
changes  in  latency  that  would  invalidate  the  selection  from  one  level 
to  another.  3)  The  predictive  power  of  a DF  obtained  from  data  of 
higher  intensity  and  applied  on  data  from  lower  intensity  is  greater 
than  the  predictive  power  of  a DF  obtained  vice  versa.  This  again  ties 
up  with  the  signal-to-noise  ratio  hypothesis. 


61 


QUESTION  6;  Runs:  C-7 , C-8,  C-9,  C-10,  C-11,  C-lla 


Effectiveness  of  averaged  data  to  train  the  D.F.  as  opposed  to  raw  data. 

These  runs  were  made  in  an  effort  to  explore  the  advantages,  if  , to 
be  gained  by  using  short  averages  of  E.R.  data  to  train  a D.F. 

The  reasoning  behind  this  was  the  standard  assumption  that  averaging  will 
eliminate  irrelevant  individual  characteristics  of  single  epoch  evoked 
responses.  As  a D.F.  System  characterizes  the  aggregrate  similarities 
within  groups  and  aggregrate  differences  between  groups,  individual  fluc- 
tuations in  single  evoked  responses  may  be  reduced  to  advantage  before 
calculating  the  D.F. 

A data  set  comprising  60  epochs  each  of  Red,  Green  and  Blue  was  used  in 
runs  C-8  and  C-9.  The  averages  were  limited  to  6 epochs  each,  thus  giving 
a total  of  only  10  A.E.R.'s  (averaged  evoked  response)  per  color. 

Several  combinations  such  as  training  on  averages  and  testing  on  raw  data 
were  tried  and  the  results  are  summarized  below.  The  regularity  of  the 
averages  led  to  perfect  discrimination  of  the  training  set.  Testing  re- 
vealed a moderate  but  definite  advantaqe.  This  approach  deserves  further 
study  to  determine  an  optimal  trade-off. 


TABLE  11 


ADVANTAGES  OF  USING  A, E.R.  TO  TRAIN  A D.F. 
SYSTEM  TO  APPLY  ON  SINGLE  E.R. 


R 

G 

B 

R 

G 

B 

R 

29 

1 

0 

F 

R * 

10 

0 

ni 

G ! 

i 

j 

4 

26 

0 90.007, 

G * 

0 

10 

0)1 00.007 

B 

i 

3 

1 

26  j 

B 

0 

0 

10  ) 

R2  1 

24 

4 

2: 

r2~ 

22 

7 

1 

G2  | 

A 1 

5 

21 

4 63.337 

G2 

5 

15 

7 70.007 

B2  ! 

t 

9 

9 

12 

B2  . 

2 

2 

26 

C-8 


C-9 


(Trained  on  raw;  Applied  on  raw)  (Trained  on  averages;  ApDlied  on  raw) 
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The  discriminability  of  the  data  improved  from  90.0%  to  100.0%  and  the 
predictive  power  of  the  OF  computed  increased  from  63.3%  to  70.0%. 

The  marked  improvement  in  the  discriminability  of  the  data  as  recognized 
by  the  program  is  also  evident  in  Runs  C-7,  C-10,  C-ll,  C-lla  also  made 
on  short  averages.  Training  set  discriminability  is  perfect  in  every 
case.  It  is  reasonable  to  assume  that  the  resulting  DF  are  also  better 
performing  although  no  comparison  with  raw  DF  were  made  here. 

TABLE  12  INVESTIGATION  USING  RAW  DATA  {SINGLE  EVOKED  RESPONSES) 

TO  TRAIN  A DF  AND  ITS  APPLICATION  ON  RAW  DATA 


C-7 

C-10 

C-ll 

C-lla 

Discriminability 

of  data 

100.00 

100.00 

100.00 

100.00 

Predictive 
Power  of  DF 

43.33 

86.66 

85.00 

76.66 

QUESTION  7;  Runs:  C-l , C-2,  C-3,  C-4,  C-8,  C-12,  D-4,  D-6,  D-7,  D-8, 

D-9,  D-10,  D-l 1 , D-12,  D-21 , D-22,  D-23 


D.F.  trained  on  single  E.R.  and  applied  on  single  E.R. 


TABLE  13  PERFORMANCE  OF  D.F.  OBTAINED  WITH  RAW  (SINGLE  E.R.) 

AND  APPLIED  ON  RAW 


Run  Code 


Discriminabil  it.y  of  Data  Predictive  Power  of  DF 


The  runs  listed  in  Table  13  were  made  for  evaluating  several  different 
hypotheses;  however,  they  are  all  examples  of  training  a D.F.  system  on 
single  E.R.'s  followed  by  application  on  single  E.R. ' s 

The  discriminability  of  the  data,  as  shown  ranges  from  66.25%  to  100.00%, 
whereas  the  predictive  power  of  D.F.  varies  from  40.58%  to  92.00%.  Run  D-4 
was  found  questionable  and  discarded.  The  percentages  indicated  here  are 
not  all  comparable  as  they  refer  to  experiments  with  different  input  sets. 
(They  have  since  been  replaced  by  a standardized  measure  of  mutual  infor- 
mation). The  technique  is  clearly  quite  worthy  of  consideration.  Yet, 
training  a D.F.  on  short  averaged  E.R.  could  prove  to  be  more  efficient. 


QUESTION  8;  Runs:  B-3,  B-4,  B-5,  B-6,  B-7 


Investigation  of  variation  in  E.R.  as  the  whole  spectrum  is  spanned  with 
five  occipital  channels. 


TABLE  14  DISCRIMINABILITY  OF  E.R.  FOR  WHOLE  SPECTRUM  AS  RECORDED  BY 
DIFFERENT  CHANNELS 


% Discriminability  of  Data  (5  Channels) 


ID  Color Wavelength(nm)  Top  Central  Bottom  Left  Right 


1. 

Violet 

420 

69 

70 

72 

63 

61 

2. 

Blue 

465 

42 

53 

TT 

17 

28 

3. 

Aqua 

490 

41 

5? 

38 

38 

17 

4. 

Green 

515 

55 

re 

6 

41 

34 

5. 

Green  Yellow 

530 

re 

42 

25 

20 

48 

6. 

Yellow  Green 

550 

Yf 

13 

20 

13 

19 

7. 

Yellow 

575 

M 

50 

31 

27 

13 

8. 

9. 

Orange 
Red  Orange 

595 

620 

25 

45 

re 

re 

8 

41 

31 

52 

16 

48 

10. 

Red 

660 

77 

73 

45 

59 

50 

Overall  discriminability  recorded  by 
each  channel  {%) 

47.66 

50.16 

31.72 

35.94 

34.06 

As  we  note  in  Table  14  the  two  channels.  Top  and  Central  Occipital/Ears 
saw  the  most  discriminability  in  the  data.  We  also  find  that  Yellow-Green 
(#6)  was  the  least  discriminable  color.  Most  of  the  epochs  recorded  for 
that  group  were  confused  to  be  from  all  the  rest  of  the  group. 

The  underlined  percentages  in  Table  14  show  that  using  the  whole  spectrum 
data,  the  extremes  are  the  most  discernable. 


The  channel  which  recorded  the  most  discriminability  for  violet  was  Bottom 
Occipital/Ears,  whereas  Red  was  most  clearly  seen  by  Top  Occipital/Ears. 

We  should  note  that  the  closeness  of  the  groups,  especially  groups  3 to  7 and 
the  large  amounts  of  data  (64  ER's  per  group)  had  a large  effect  in  the  final 
discriminability  of  data. 
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QUESTION  9;  Runs:  C-3,  C-4 


To  study  the  closeness  of  Green,  Green-Yellow  and  Yellow-Green  In  comparison 
to  the  extremes  of  the  spectrum. 


The  runs  B-3,  B-4,  B-5,  B-6,  B-7  showed  us  that  responses  recorded  for 
groups  4,5,6  were  very  much  similar  (also  see  Question  8 ).  This 

prompted  a run  to  study  these  groups  in  comparison  to  the  extremes  of  the 
spectrum. 

From  this  run  we  see  that  on  the  application  of  DF  developed  on  the  above 
groups,  4,  5,  and  6 do  show  higher  discriminabillty  when  considered  by  them- 
selves; however,  group  6 is  more  like  group  4. 


QUESTION  10;  Runs:  B-3,  B-4;  B-5,  B-6,  B-7,  B-8,  B-9 
C-l,  C-2,  C-3,  C-4,  C-5 


To  investigate  wrap-around  effects  of  the  spectrum. 

Looking  at  results  obtained  from  runs  made  on  the  whole  spectrum  led 
to  speculation  that  the  effect  of  variation  in  the  wavelength  of  the  stimulus 
resulted  in  a corresponding  gradual  variation  in  ER  and  that  the  extremes 
of  the  spectrum  resulted  in  similar  ER 1 s suggesting  a wrap-around  effect. 

As  we  notice  in  Runs  B-3,  B-4,  B-5,  B-6  a large  number  of  ER's  at  one  end 
were  classified  as  coming  from  the  other  end.  The  following  are  some  such 
examples  in  Table  15, 
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TABLE  15 


SUGGESTIONS  OF  A WRAP-AROUND  EFFECT 


Correctly  Discriminated  Incorrectly 

Run  ID  Group  Classified^) as  Group  Discriminated  (%) 


B-3 

2 

17.19 

10 

15.63 

B-4 

8 

48.44 

2 

10.94 

B-5 

1 

68.75 

9 

10.94 

B-5 

8 

25.00 

2 

10.94 

B-5 

9 

45.31 

1 

10.94 

B-5 

9 

45.31 

1 

10.94 

B-6 

8 

07.81 

1 

n .88 

B-7 

1 

60.94 

10 

12.50 

B-7 

8 

15.63 

2 

15.63 

B-7 

10 

51.56 

1 

15.63 

The  incorrect  classification  numbers  are  not  extremely  larqe,  but  they  are 
the  next  highest  percentages  for  those  Individual  groups  and  thus  gain  sig- 
nificance. As  a result,  runs  C-l,  C-2,  C-3,  C-4  and  C-5  were  made  with 
data  belonging  to  the  extremes  to  see  if  the  DF  system  trained  on  such  data 
would  incorrectly  classify  a significant  amount  of  ER's  by  putting  them  Into 
opposite  camps.  This  did  not  occur  and  on  the  contrary  the  extreme  groups 
were  found  to  be  well  apart  from  each  o^her. 
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QUESTION  11;  Runs:  B-3.  B-4,  B-S,  B-6,  B-7,  B-8,  B-9 


To  study  the  performance  of  DF  trained  on  Red,  Green  and  Blue  ER's,  and 
applied  on  ER's  recorded  due  to  the  rest  of  the  spectrum. 


These  runs  were  made  to  Investigate  the  possibilities  of  being  able  to 
describe  ER's  recorded  by  the  spectrum  in  terms  of  ER's  recorded  by  Red, 
Green  and  Blue. 

This  was  accomplished  by  training  a DF  system  on  Red,  Green  and  Blue,  and 
studying  its  performance  in  predicting  ER's  belonging  to  the  rest  of  the 
spectrum. 

By  looking  at  the  final  classification  matrices  of  Runs  B-8,  B-9  no  obvious 
trends  were  found. 


QUESTION  12;  Runs:  D-11,  0-12,  D-22 


Investiqation  of  the  discriminability  of  E.R.'s  obtained  with  opponent 
color  pairs. 


The  following  statistics  were  recorded  when  runs  of  opponent  color  pairs  were 
made.  The  opponent  color  pairs  were  chosen  by  usinq  the  chromaticity 
diaqram  for  standard  source  "C"  (artificial  dayliqht).  (See  Figure  1). 

The  opponent  color  for  each  pair  was  selected  by  extendinq  the  line  con- 
nectina  the  color  and  White  to  the  opposite  side. 


Opponent 
Color  Pair 

Discriminability 
of  Data  [%) 

Predictive  Power 
of  DF  (%) 

Red-Blue-Green 

100.00 

88.00 

Green-Purple 

98.00 

78.00 

Blue-Yellow 

98.00 

8b.  00 

We  note  that  while  the  discriminability  recorded  by  the  Stepwise  Discriminant 
Analysis  proqram  was  excellent  for  all  these  opponent  color  pairs,  the  per- 
formance of  the  DF  system  computed  for  the  Green-Purple  pair  was  siqnificantly 
lower  but  still  on  the  hiqh  side  when  compared  to  typical  runs  in  other  com- 
binations . 

Examininq  the  sample  points  selected  in  these  runs  did  not  lead  at  this  time 
to  any  conclusive  arguments  in  favor  of  component  codes  of  the  Red,  Green 
and  Blue  processes. 


QUESTION  13;  Runs:  D-19,  D-21 


To  study  the  discriminabil 1 ty  of  opponent  colors  of  Red,  Green  and  Blue, 
i.e.  Blue-Green,  Purple  and  Yellow  respectively. 


In  conjunction  with  Question  12  this  run  was  made  to  study  the  dis- 
criminabil ity  of  Blue-Green,  Purple  and  Yellow.  The  results  were  very 
much  similar  to  those  obtained  for  Red,  Green  and  Blue.  They  are  as 
fol lows : 


Discriminabil  ity 

Run  ID  Color  Tripole of  Data  (%) 


Predictive 
Power  of  DF  {%) 


D-19  BG,P,Y 

D-21  R,  G,  B 


90.67  73.33 

91.67  77.08 
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QUESTION  14:  Runs:  C-13,  C-13a,  C-14,  C-15,  C-16,  C-18,  C-19 


Search  for  component  codes  by  discriminating  against  a constant  DC  response. 

The  rationale  behind  this  experiment  was  to  compare  the  responses  with 
respect  to  a neutral  "group".  The  results  hopefully  would  identify  the 
specific  time  windows  where  the  activity  from  individual  colors  was  at 
a peak. 

The  following  table  summarizes  the  resulting  windows  In  the  order  deter- 
mined by  the  Stepwise  Discriminant  Analysis  program. 


TABLE  16  TIME  VARIABLE  DISTINCTIVES  OF  RED,  GREEN,  BLUE, 

YELLOW,  ORANGE  & VIOLET 


Run  Code Color  Variables 


C-l  3a 

Red 

14, 

20, 

22, 

28, 

44, 

25, 

55, 

63, 

18, 

15 

C-14 

Green 

14, 

17, 

38, 

34, 

43, 

F0, 

29, 

31, 

16, 

64 

C-16 

Blue 

53, 

16, 

63, 

21. 

13, 

26, 

27, 

58, 

79, 

69 

C-l  3 

Yel low 

25, 

62, 

16, 

43, 

60, 

48, 

34, 

41, 

37, 

44 

C-19 

Violet 

13, 

35, 

18, 

46, 

32, 

11, 

73, 

14, 

48, 

65 

C-18 

Orange 

35, 

59, 

12, 

42, 

64, 

32, 

13, 

61, 

20, 

14 

The  adjoining  Figure  is  a plot  of  the  above  variables  in  relative  order 
for  the  Red,  Green,  Blue  and  Yellow  colors.  A wrap-around  effect  from  Red 
to  Blue  is  visible.  It  must  be  noted  that  only  time  points  are  considered, 
regardless  of  amplitude. 


",  p»r 

\\  \h  Si 


r.j 


\ .J  t 

m - 


^•PLDT:>e  sx  sy  -?: 


y 

•» 1< . 


■i  '?r*Tr  IS  THE  T1TUE?  BCI  COLOR-. VERSUS  CONSTATN««NT 

It  THE  L**EL  FDR  the:^rxzs?  sample  point  number 
, Whrt  is  the  label  for  the  yaxis?  coo^lor  code 
v What  is  the  symbol?  ♦ 

>*#Wmat.  is  the  lower  limit  ■ 


• ' V w*‘  ’ 

$ r: 

.1- 


K^>C*  FOR  THE  XT  axis?  1 
%What,  is  the.  upper  limit  for  THE  X AXIS?  80 


MM&;  * r^jg 


AXIS?  1 


.-  » { • - \ * V?*  v 


lS  THE  lower  LIMIT  FOR.  THE  Y 
‘(^What  is  the  upper  limit  for. the 
'”>N'CY-~  AXIS?  5 

Whatms  the  page  width  iil characters?  (Must  be  at  least  ?£>  75 
|^|“,ORUND>NNnTHE’‘ S"ME  S“LjDN  *°™  **«*?••  <VDL,  HUST  KNSHER  VE^ 
Bcr^  COLOR  versus  constant 


..." 

Vis  •sX'Sfe'V4*-*. 

> &£%&“’  - . . fdsW’KJSW* 


W:  ■ 


rf  ♦*«'  COLOR  CODE 

J'^5.0000 


. f-  w “V  «r  ..... 

.;•  ....  • 

• '■si : : 

. -,r*  . 


T 


4. 00  00 


.$£$77  ‘ 


&.0000 


aXoooo  + 


V '•♦ — 


.T 

V 


I * 


v'  ::} 


♦ ♦ ♦ ♦ ♦ 


♦ ♦ 


YELLOW 


♦ ♦-  ♦ 

s'. tj..'''rfi 


■-  I 

*‘  • V*.‘-  VrV 
-,;..Ain;,r4*rir  -.• 

• /•**'».'<'_  I* 

■•  ■ i .»*'  C ' ** 


♦ • ♦- 


• M. 


I' 


.BLU0  . ..... 

%;  ....-35 

V J : 4 * *,  •'  - V ~-Zi^ 

Ttfsfcs.:: 

^ .-GREEN  •*>;  MJ:wV 


>.0000 

■V: 


•Tv.i 




0.0000 


♦ ♦-  ♦-  ♦- 

+. 


£0.0000 


-+- 

40.0000 


-+-- 

60.0000 


+-  :•  .v  ’ 


80.0000  'v  %%  -3 


T R^l.564 


i . 

. 4E1g. 


25 


SAMPLE  POINT  NUMBER 

* l»  •,.  * * ■'• 

Discrimination  Aqainst-WHITE  ER'l 


-t*/: 


.Tfvvr'-i  K '■*  • *■  .'1.1 


•*vv,r- 


■> ; j;-T 

i J*‘«. 


a 


72 


■ 


0 


QUESTION  15;  Runs:  D-6,  D-7,  D-8,  D-9,  D-10 


In  search  for  component  codes  by  discrimination  against  ER's  recorded 
for  White  color. 


Following  up  some  of  the  search  strategy  as  in  Question  14,  these  runs 
were  made  to  compare  ER's  recorded  for  Red,  Green,  Blue,  Yellow  and  Purple 
discriminated  against  a group  of  ER's  recorded  with  a White  flash. 


TABLE  17  TIME  VARIABLES  CHOSEN  FOR  DISCRIMINATION  AGAINST 

WHITE  COLOR 


Run  Code 


Color 


Variables  Chosen 
(in  order) 


v- 


D-10 

Red 

32, 

9, 

42, 

49, 

72, 

68 

D-9 

Green 

66, 

36, 

57, 

26, 

68, 

64 

D-8 

Blue 

54, 

64, 

22, 

61, 

33, 

75 

D-7 

Yellow 

63, 

65, 

72, 

70, 

56, 

38 

D-6 

Purple 

25, 

75, 

55, 

24, 

63, 

72 

The  discriminability  of  the  data  was  again  excellent.  The  predictive 
power  of  DF  was  lower  in  the  Yellow  run  which  was  against  a Yellow 
background,  and  is  also  less  discernable  behavioral ly. 


A chart  similar  to  Fig.  25  is  shown  in  Fig. 26  where  variables  chosen 
for  each  color  against  White  are  plotted  along  with  a subset  of  variables 
chosen  against  a constant,  pointing  to  a fairly  complex  relation- 


nst  WHITE 
nst  Constan 


Responses 


QUESTION  16;  Runs:  C-10,  C-11 , C-lla 


To  investigate  the  possibility  of  expressing  ER's  due  to  one  color  in 
terms  of  the  other  two  in  the  case  of  tripole  of  Red,  Green  and  Blue. 

An  attempt  was  made  to  study  the  similarities  between  ER's  due  to  the 
three  colors  namely.  Red,  Green  and  Blue.  In  these  runs  the  DF  was 
trained  on  the  AER's  from  two  groups  and  applied  on  AER  and  single  ER 
from  the  third  group.  The  results  are  as  follows: 


lun  Cod< 

Training  Testing 

2 Groups  Group 

Classified  (%)  as 

Group  1 Group  2 

All  Raw  All  Raw 

Final 

F-value 

C-10 

1.  Red  2.  Blue 

Green 

80 

63 

20 

37 

914.23 

C-11 

1.  Green  2.  Blue 

Red 

30 

27 

70 

73 

1095.90 

C-lla 

1 . Red  2.  Green 

Blue 

10 

40 

90 

60 

1450.26 

An  interesting  circulating  symmetry  is  observed  as  shown  in  Fig. 


Fig.  27  DF  Trained  on  Two  Groups  of  Tripole 
and  Applied  on  the  Remaining  Third 


Red 


In  Fig.  27  the  arrow  Indicates  the  direction  of  the  classification  of 
of  each  color,  given  a choice  between  the  two  others  used  to  create  the 

DF  as  the  table  above  shows.  Note  the  unusually  large  final  approximate 
F-values. 


QUESTION  17:  RurK  D-23 


To  evaluate  the  window  estimated  by  runs  associated  with  Red,  Green  and 
Blue  colors. 

One  of  the  most  important  results  obtained  In  this  phase  of  data  analysis 
has  been  the  Identification  of  the  time  windows  recognized  by  Stepwise 
Discriminant  Analysis  as  containing  the  most  information  reflecting  the 
separability  of  ER's  recorded  for  the  Red,  Green  and  Blue  stimuli. 

Table  18  sunmarlzes  the  time  points  selected  (In  msec.)  in  the  more 
significant  runs.  A simple  histogram  (Fig. 28  ) of  these  time  points 
shows  that  most  of  the  Information  seems  contained  in  time  windows  aroun 
88,  100,  116,  132  and  212  msec,  taken  with  an  aperture  of  + 4 msec. 

This  result  is  important  because  It  provides  us  a basis  for  dimensionality 
reduction  by  discarding  the  other  time  points  from  the  analysis  in 
color  discrimination.  The  data  reduction  is  a precious  asset  to  study 
DF's  trained  on  multi-channel  data  because  of  the  reduction  in  the  number 

of  variables  to  consider. 


An  evaluator*  discriminant  analysis  run  was  made  using  the  variable  windows 
of  23-29,  34-38,  43-48  and  51-56. 

As  Run  D-23  shows,  the  dlscrlmlnabllity  of  data  was  maintained  at  78.6%, 
and  at  64%  the  predictive  power  of  DF  did  not  suffer  any  appreciable  loss. 
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TABLE  18  GENERAL  AREAS  OF  ACTIVITY  FOR  RED,  GREEN  AND  BLUE  STIMULI 


Run  Code  Channel  § Exp.  ID. 


Variables  Selected  (First  Six 


9 

1 

5 

2 

3 

Three,  Color 


Ten  Color 


Opponent  Color 


50,  69, 
50,  40, 
41,  34, 
20,  44, 
53,  44, 

52,  73, 
74,  37, 
39,  37, 
74,  37, 
73,  57, 

53,  44, 
53,  44, 
53,  44, 
53,  44, 
15,  63, 

52,  20, 

53,  18, 
53,  18, 
27,  53, 
30,  66, 
36,  25, 
53,  18, 
53,  31, 
25,  33, 


61,  41, 
69,  61, 
45,  29, 
13,  9, 

62,  29, 
54,  35, 
39,  52, 

58,  75, 
39,  54, 
51,  9, 
62,  29, 
74,  48, 
74,  48, 
62,  29, 
53,  42, 
76,  28, 
25,  36, 
25,  36, 
66,  33, 

59,  22, 
53,  32, 
25,  36, 
36,  18, 
66,  47, 


32,  30 
36,  30 
26,  22 

32,  26 
67,  48 
38,  63 
18,  72 
48,  57 
72,  68 
53,  33 
67,  48 
67,  62 
42,  38 
67,  48 
21,  59 
62,  25 

33,  22 
33,  22 
70,  24 
33,  36 
48,  69 
33,  22 

25,  50 

26,  29 


QUESTION  18:  Runs:  D-l , D-2»  D-'2 

To  investigate  the  effect  on  the  predictive  power  of  DF  of  the  a priori 
elimination  of  statistical  outliers  from  the  computation  of  DF. 

. ser,es  of  three  Iterative  runs  was  made  to  see  if  elimination  of 
"outliers11  (responses  which  ware  substantially  outside  the  cluster  of 
their  qroup)  enhanced  the  discrlminability  of  data. 

Table  19  below  sumarizes  the  results  obtained  from  these  runs. 

ITERATIVE  REJECTION  OF  OUTLIERS  

— Distances  between  means  of 

Discrlminability  Final  clusters  on  th< i canonical 
Run  Code  Iteration#  of  data  («  U-.Utirtic^  Red-B.lue 


91.67 

95.45 

100.00 


0.084 
0.C  8 
0.026 


5.229 

5.987 

6.402 


2.641 

4.201 

4.692 


4.25 

4.902 

5.613 


,t  can  oe  seen  that  the  discrlminability  of  data  has  improved  con  id.  ably 
with  each  iteration.  The  results  are  Illustrated  in  the  canonical  map 
Of  the  three  runs.  The  clusters  of  Red.  Green  and  Blue  not  only  trave 
farther  from  each  other,  but  also  be.ome  more  condensed. 

The  iterative  elimination  of  "outliers"  from  the  training  set  is  now  con- 
sidered in  the  methodology.  It  must  be  noted  that  the  variables  selec 
in  the  three  runs  were  almost  identical  thus  giving  further  credence 
the  validity  of  the  DF  system  computed.  


Run  Code 


Variables  Selected 

53,  44,  62,  29,  67,  48,  38,  42,  75,  51 

53,  44,  74,  48,  67,  62,  38,  42,  51,  33 

53,  44,  74,  48,  42,  38,  67,  41,  62,  33 
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SIMULATION  RUNS  OF  THE  ORTHOGONAL  DERIVATION 
PROGRAM 

Tests  Results  on  Ortho 

The  orthogonal  derivation  program  (ORTHO)  was  subjected  to  a battery 
of  tests  on  the  synthetic  data.  Two  new  programs  were  written  for  that 
purpose.  One  is  a general  purpose  PRINT-PLOT  routine  and  the  other  one 
is  a synthetic  data  generator  (SDG) . SDG  was  used  to  qenerate  three  sets 
of  data  each  consisting  of  3 channels  with  1360  samples. 

Data  Set  1 (f  = 4Hz) 

Xj  = random  nois?  with  uniform  distribution  (RN) 
x2  = Xj  + sin  (2rr*4*t) 
x3  = x2  for  1 , 1280 

x2  + VER  for  1281  , 1360 

where  VER  = s'm  (2tf*l 0*t) 

Data  Set  2 (f  = 4Hz) 

Xj  = RN 

x2  = x^sin  (2tf*4*t) 
x3  = x2  for  1 , 1280 

x2  + VER  for  1281,1360 

where  VER  is  the  same  as  in  Data  Set  1 
Data  Set  3 (f  = 10Hz) 

The  same  as  the  Data  Set  1,  except  f = 10Hz 

The  tests  performed  were  intended  to  answer  the  followinq  questions: 

c.)  will  ORTHO  separate  VER  from  data  set  1 and  2?(i.e.  addition 
and  multiplication) 

b)  how  ORTHO  behaves  with  changing  frequency  (4hz  and  10Hz  were 
used)? 
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c)  how  many  samples  (time)  are  needed  for  VER  to  qet  a qood  separa- 
tion from  VER? 


The 

following  tests  were  performed: 

ORTHO  was  Run  on 

Data  Set  (order  of 

is  always 

1. 

full  length  of  epoch 

(f  = 4Hz,  DS  1 , Xj  x.  x3) 

2. 

full  length  of  epoch 

(f  = 4Hz , DS  2,  Xj  x2  x3) 

3. 

full  length  of  epoch 

(f  =10Hz , DS  3,  Xj  x2  x3) 

4. 

full  length  of  epoch 

(f  = 4Hz , DS  1,  x3  x2  Xj) 

5. 

70  sample*  before  VER 

(f  = 4Hz,  DS  1 , Xj  x2  x3) 

6. 

70  samples  before  VER 

(f  =10Hz , DS  3,  x1  x2  x3) 

7. 

130  samples  before  VER 

(f  = 4Hz , DS  1 , x-j  x2  x3) 

8. 

130  samples  before  VER 

(f  =10Hz , DS  3,  x1  x2  x3) 

Preliminary  Conclusions 

ORTHO  separated  VER  very  well  on  all  data  sets,  only  one  difference 
being  that  when  full  length  of  epoch  was  used  the  amplitude  of  separated 
VER  was  about  twice  as  large  as  if  only  130  or  70  samples  were  used. 

Thus  Is  appears  that  with  the  present  filter  constants  only  about  100 
samples  before  VER  are  needed  to  qet  good  separation  of  VER. 
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white  noise 


white  noise  + 10 hz  sine  umuc 


STEPWISE  DISCRIMINANT  ANALYSIS  RUN 
SUMMARIES 
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BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH4RAW 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

SAMPLING  RATE  - 4 MSEC. 


IHANNEL  PROCESSED:  4-Eye  (Inferior  Orbital)  siJBJ EC T'S  INITIAL:  XXX 

Ears 

J UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Red,  48/Green,  47 /Blue  tested) 


F-LEVEL  FOR  INCLUSION:  0.0 


O.  OF  STEPS  RUN:  82 


F-LEVEL  FOR  DELETION:  0.0 


MEMO: 


The  following  results  are  after  10  steps 


RESULTS  OF  RUN:  A-l 


VARIABLES  SELECTED:  (in  order) 

20,  44,  13,  9,  32,  26,  29,  4,  12,  48 


>ERFORMANCE 

Mscriminability 

of  Data 

88  . 81  o/o 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
-STAT  0.54687 


FINAL 

0.11940 


Approx.-  F 58.001341  24.81071 


^INAL  CLASSIFICATION  MATRIX: 

I R G 


R 

45 

1 

2 

G 

1 

41 

6 

B 

1 

5 

41 

85 

Yellow 


Blue-G 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH3NEW 


NO.  OF  TIME  POINTS:  80 

HANNEL  PROCESSED:  3-Occipital/Ears 


DATA  CODE: 


SAMPLING  RATE  4 MSEC, 


HANNEL  PROCESSED:  3-Occipita.l/Ears  SlJHJECTfS  INITIAL*  XXX 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  (•-'  indicates  group  is  being 
24/Red,  24/Green,  24/Blue,  tested) 

-24/New  Red,  -24/New  Green,  -23/New  Blue 

F-LEVEL  FOR  INCLUSION:  0.01  F-LEVEL  FOR  DELETION:  0. 005 


>0.  OF  SI  EPS  RUN:  10 


MEMO: 


RESULTS  OF  RUN:  A-2 


VARIABLES  SELECTED:  (in  order) 

32,  41,  31,  47,  17,  12,  44,  19,  22,  25 


PERFORMANCE 
Mscriminability  Predictive 

of  Data Power  of  DF 

100*  00  % 85  • 92  o/o 


INITIAL 

J-STAT  0.51348 


FINAL 


0.04684 


Approx.-  F32.68846  21.72443 


INAL  CLASSIFICATION  MATRIX: 


R 

24 

0 

0 

G 

0 

24 

0 

B 

0 

0 

24 

New  Red 

18 

0 

6 

New  Green 

0 

24 

0 

New  Blue 

4 

0 

19 

86 

Y ellow 


Blue-G 


_ i 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH3RAW 


"0.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAM PL1NG  RATE  4 M SEC . 


HANNEL  PROCESSED:  3- Occipital/Ears 


SUBJECT'S  INITIAL:  XXX 


UMBER  of  Cases/Names  of  Groups  Analyzed:  ('-'  indicates  group  is  being 
48/Red,  48/Green,  47/Blue  tested) 


[F-LEVEL  FOR  INCLUSION:  0.0 
O,  OF  STEPS  RUN:  10 


F-LEVEL  FOR  DELETION:  0.0 


MEMO: 


RESULTS  OF  RUN: 


MRIABLES  SELECTED;  (in  order) 

41,  34,  45,  29,  26,  22,  33,  10,  68,  65 


PERFORMANCE 
Jiscrim  inability 
of  Data 
92  . 31°/o 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
J-STAT  0.53200 


FINAL 

0.06551 


Approx.-  f 61. 578861  38.08063 


INAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

R 

44 

0 

4 

G 

0 

48 

0 

B 

7 

0 

40 

Red 

Grcjn 

Blue 

Yellow 

White 

Purple 

Blue-Gi 


BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  BCI07M 


O.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAM  PUN  G RATE  4MSEC. 


:HANNEL  PROCESSED:  3-Occipital/Ears  SUBJECT'S  INITIAL:  XXX 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Red,  48/Green,  47/Blue  tested) 


-LEVEL  FOR  INCLUSION:  0.0 
O.  OF  STEPS  RUN:  25 


tested) 

F-LEVEL  FOR  DELETION:  0.0 


MEMO: 


The  results  are  after  25  steps 


RESULTS  OF  RUN:  A-4 


VARIABLES  SELECTED:  (in  order) 

41,  34,  45,  29,  26,  22,  33,  10,  68,  65, 
36,  71,  79,  24,  21,  17,  20,  25,  27,  16, 
70,  73,  39,  38,  80 


>ERFORMANCE 

Mscriminability 

of  Data 

lOOpO  o/o 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
-STAT  0.53200 


FINAL 


0.02622 


Approx.-  F 61. 57886  24.01642 


•INAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

R 

48 

0 

0 

G 

0 

48 

0 

B 

0 

0 

47 

88 

Yellow 


Blue-G 
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BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  C3GRBL 


NO.  OF  TIM K POINTS:  80 


DATA  COOK: 


SAMPLING  HATK  1MSEC 


HANNEL  PROCESSED:  3-Occipital/Ears 


SUBJECT’S  INITIAL:  XXX 


1UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Green,  47/Blue  tested) 


F-LEVEL  For'  INCLUSION:  0.0 

O.  OF  STEP.  .\ UN:  10 


F-LEVEL  FOR  DELETION:  0.0 


MEMO: 


RESULTS  OF  RUN;  j 


VARIABLES  SELECTED:  (in  order) 

41,  34,  45,  10,  37,  64,  73,  80,  59,  18 


ERFORMANCE 


[scrim  Inability 
of  Data 


100  • 00  °/° 


Predictive 
Power  of  DF 


INITIAL 


J-STAT  0-40792 


FINAL 


Approx,-  F 134.  9834 


INAL  CLASSIFICATION  MATRIX: 


G 

48 

0 

B 

0 

47 

Key 

R 

Red 

C = 

Green 

B = 

Blue 

Y = 

Y ellow 

W = 

White 

P = 

Purple 

89 

BG  = 

Blue-Gi 

I 


r 


I 


f 


BMP  07  M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  C3  REGR 


JO.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


IHANNEL  PROCESSED:  3 -Occipital/Ears 


SUBJECT'S  INITIAL:  XXX 


1 UMBER  of  Cases/Names  of  Groups  Analyzed;  indicates  group  is  being 
48/Red,  48/Green  tested) 


F-LEVEL  FOR  INCLUSION:  0.0 
O.  OF  STEPS  RUN:  10 


F-LEVEL  FOR  DELETION:  0.0 


RESULTS  OF  RUN;  A-8 

VARIABLES  SELECTED:  (in  order) 

33,  40,  10,  46,  29,  41,  21,  73,  64,  36 


PERFORMANCE 
Discrim  inability 
of  Data 

100*00  °/° 


Predictive 
Power  of  DF 
. o/o 


INITIAL  FIN  A 

U-STAT  0.57485  01682 

Approx.-  F69. 51984  64.25989 


INAL  CLASSIFICATION  MATRIX: 


Yellow 


Blue-G 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH3RED 

NO.  OF  TIME  POINTS:  80 


DATA  CODE- 

SAM  PUN C RATE  * MSEC. 


-HANNEL  PROCESSED:  30-Occipital/Ears  SUBJECT'S  INITIAL-  XXX 

i UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Red,  95 /No nr ed  tested) 


F-LEVEL  FOR  INCLUSION:  0.0 
NO.  OF  STEPS  RUN; 


F-LEVEL  FOR  DELETION:  0.0 


MEMO: 


Here  Group  'nonred'  consists  of  48/Green  and  47/Blue 


RESULTS  OF  RUN:  A-9 

VARIABLES  SELECTED:  (in  order) 

34,  29,  33,  26,  69,  38,  55,  51,  47,  36 


PERFORMANCE 
Dlscriminability 
of  Data 


93  ’ 71  °/° 


Predictive 
Power  of  DF 


INAL  CLASSIFICATION  MATRIX: 


INITIAL 


J-STAT  082891 


FINAL 


0.34836 


Approx.-  F29. 10364  24.69171 


Nonred 


R 

45 

3 

Nonred 

6 

89 

91 

BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


jNO.  OF  TIME  POINTS:  80 


CHANNEL  PROCESSED:  3-Occipital/Ear3 


RUN  CODE:  CH3GRN 


DATA  CODE; 


SAMPLING  RATE  - 4 MSEC. 


SUBJECT'S  INITIAL:  XXX 


[NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Green,  95/nongreen  tested) 


F-LEVELFOR  INCLUSION:  0.0 


F-LEVEL  FOR  DELETION: °»° 


[NO.  OF  STEPS  RUN;  10 


MEMO: 


Here  group  'Nongreen'  consists  of  48/Red  and  47 /Blue 


I 

! 


RESULTS  OF  RUN:  A-10 


[VARIABLES  SELECTED:  (in  order) 


41,  34,  45,  10,  73,  64,  37,  20,  33,  28 


[PERFORMANCE 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

100  -oo  °/° 

. o/o 



FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0 58025 

0.13568 

Approx.-  F 101.  99911 

84.08607 

G 

Nongreen 

G 

43 

0 

Nongreen 

0 

95 

92 

Key 
R = 

Red 

G = 

Green 

B = 

Blue 

Y = 

Y ellow 

N = 

White 

P = 

Purple 

BG  = 

Blue-G: 

$ 


EMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH3BLU 


DATA  COOK: 


MEMO: 


1.  Here  the  group  'Nonblue'  consists  of  48 /Red  and  48/Green 

2.  The  following  results  are  after  10  steps. 


RESULTS  OF  RUN:A-11 

fARIABlA-T.  qELECTED:  (in  order) 

40,  30,  26,  67,  33,  22,  38,  24,  21, 

Note:  On  step  10,  variable  26  was  removed  as  it  had  a F-value  of  0.2100 


PERFORMANCE 

Mscriminability 

of  Data 

95  . 11  o/o 


Predictive 
Power  of  DF 
. o/o 


INITIAL 

U-STAT  0.72677 


FINAL 

0.35493 


Approx,-  F53.00856  30.44206 


•’INAL  CLASSIFICATION  MATRIX; 


Nonblue 


Nonblue 


Key 

Ik  = Red 

m = Green 

IF  = Blue 

V Yellow 
W = White 

P = Purple  I 

BG  = Blue-Grn. 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  C3REBL  DATA  C()DK. 

NO.  OK  TIM E POINTS:  80  SAMPLING  RATE  4 MSEC. 

CHANNEL  PROCESSED:  3-Occipital/Ears  SUBJECT'S  INITIAL:  XXX 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

48/Red,  47/Blue  tested) 


F-LEVEL  FOR  INCLUSION:  0.0 
NO.  OF  STEPS  RUN;  10 

MEMO: 


F-LEVEL  FOR  DELETION:  0.0 


RESULTS  OF  RUN:  A-12a 

VARIABLES  SELECTED;  (in  order) 

29,  33,  26,  22,  79,  24,  21,  17,  20,  16 


PERFORMANCE 

Discriminability 

of  Data 

°/° 


Predictive 
Power  of  DF 
. o/o 


97  90  [ , ’ 

FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.81549 

0.193&7 

Approx,-  F21-04158 

32. 92717 

B 

0 

Key 

R = Red 
G = Green 
B = Blue 
Y = Y ellow 
W = White 
P = Purple 
BG  = Blue-Grn. 


BMP  07  M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CH2RAW 


0 


DATA  CODE: 


NO.  OF  TIME  POINTS;  80 

SAMPLING  RATE  4 MSEC. 

CHANNEL  PROCESSED;  2-Vertex/Eare 

SI!  UJ  LOT'S  INITIAL  - XXX 

NUMBER  of  Cases/Names  of  Groups  Analyzed: 

indicates  group  is  being 

tested) 

48  Red,  48 /Green,  47 /Blue 

F-LEVEL  FOR  INCLUSION:  n n 

F-LEVEL  FOR  DELETION:  n 0 

NO.  OF  STEPS  RUN:  81 

MEMO: 


The  following  results  are  after  15  steps 


0 

RESULTS  OF  RUN:  A-14 

VARIABLES  SELECTED:  (in  ordei 

•) 

50,  40 

, 69,  61,  36,  30,  44,  74,  65, 

80,  42,  38,  56,  32, 

27 

0 

PERFORMANCE 

INITIAL 

FINAL 

0 

Discrimina^  ility 

Predictive 

U-STAT  0.63463 

0.19242 

of  Data 

Power  of  DF 

. o/o 

. o/o 

^ 

Approx.-  F 40.2  9965 

10.74932 

FINAL  CLASSIFICATION  MATRIX: 

— i 

0 

R 

G B 

44 

' 

n 

2 2 

Key 

G 

1 

41  6 

R 

= Red 

B 

G 

= Green 

2 

4 41 

3 

= Blue 

Y 

= Y ellow 

W 

= White 

P 

- Purple 

• 

95 

BG 

= Blue-Gi 

f 

BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  COPE:  CH2RAW 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

SAMPLING  RATE  -4  ML.EC. 


H ANN  LI  PR  )C  hSSED;  2— Verten/Ears  SUBJECT'S  INITIAL"  XXX 

J UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 


48 /Red,  48,  Green,  47/Blue 

F-LEVEL  FOR  INCLUSION:  0.0 
O.  OF  STEPS  RUN:  10 


MEMO: 


tested) 

F-LEVEL  FOR  DELETION:  0.0 


RESULTS  OF  RUN;  A-15 
VARIABLES  SELECTED:  (in  order) 

50,  40,  69,  61,  36,  30,  44,  74,  65,  80 


PERFORMANCE 
Hscriminability 
of  Data 

83  • 92  °/° 


Predictive 
Power  of  DF 
. o/o 


INITIAL FINAL 

U-STAT  0.63463  0.24954 

Approx.-  F 40.29965  13.12502 


FINAL  CLASSIFICATION  MATRIX: 

‘ _J 

R G B 

R 

4|  3 4 

G 

3 41  4 

B 

4 7 36 

56 

DMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CII1RAW 


O.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


HANNE1.  PROCESSED:  1-Frontal/Ears 


SUBJECT'S  INITIAL:  XXX 


JUMPER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

. . tested) 

48/Red,  48/Green,  47/Blue 


F-LEVEL  FOR  INCLUSION:  0.0 


O.  OF  STEPS  RUN:  10 


MEMO: 


F-LEVEL  FOR  DELETION:  0.0 


RESULTS  OF  RUN:  A-18 


/ARIABLES  SELECTED:  (in  order) 

50,  69,  61,  41,  32,  30,  20,  36,  4,  40 


ERFORMANCE 


Discriminability 
of  Data 

Predictive 
Power  of  DF 

. o/o 

. o/o 

-J 

INAL  CLASSIFICATION  MATRIX: 


INITIAL 
■3TAT0.  77242 


FINAL 


0.28146 


Approx.-  F20. 62447  11.59229 


Y ellow 


Blue-G 


I 

1 
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BMD07M  STEPWISE  DISCRTMWakt 


RUN  CODE;  TESTRN 
{NO.  OF  TIM E POINTS:  80 


DATA  CODE: 

SAM  PUNK  HATri  4 MSEC 


CHANNEL  PROCESSED:  CHI-Top  Occipital/Ears  „ 

CHZ -Central  OcctpiM/Ears  SU1!JKC,S  Wrr“L:  CP 

NWIBER  of  Cases/Names  of  Groups  An.iv,„a.  indlcatcs  group  ls  Mng 

32/1,  32/5,  32/A,  tested) 

-32/U,  -32/F,  -32/T 

F-LEVEL  FOR  INCT.IISTnu.  0 010  - T mr, 

' * U F-LEVEL  FOR  DELETION- 0. 005 

NO.  OF  STEPS  RUN. 


MEMO: 


1.  Only  the  first  10  variables,  i.  e.  first  40  MSU 
were  considered  in  this  analysis. 

2.  The  odd  variables  were  taken  from  CHI,  and  the 
even  variables  were  taken  from  CH2. 


RESULTS  OF  RUN:  B-l 
VARIABLES  SELECTED:  (in  order) 

Channel  1 


1.  3.  8,  7,  9,  4,  2,  10,  5 
Channel  2 


performance  — 

Discriminability  Predictive 

of  Data Power  of  DF 

47  • 92 0/0  34.38  o/o 

FINAL  CLASSIFICATION  MATH fX- 


INITIAL 

J-STAT  0.96796 


final 

0.84322 


Approx,-  f 1.53907  I 0.74766 


1 

5 

A 

15 

8 

9 

9 

15 

8 

11 

5 

16 

10 

11 

11 

10 

10 

12 

11 

8 

13 

[ l=Violet 
5=Green- 
Yellow 
A=Red 
U=Violet 
5=Green- 
Yellow 
T=Red 


Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-Grn. 


c 


0 


o 


c 


K 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


(no.  OF  TIME  POINTS:  80 


CHANNEL  PROCESSED;  4-Left  Occipital/Ears  SUBJECT'S  INITIAL:  CP 


RUN  CODE:  TESTRN 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


piMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

64/1,  64/2,  64/3,  64/4,  64/5,  tested) 


64/6,  64/7,  64/8,  64/9,  64/A 
[F-LEVEL  FOR  INCLUSION:  0.010 


F-LEVEL  FOR  DELETION:  0.005 


(NO.  OF  STEPS  RUN-  10 


MEMO: 


The  first  10  variables,  ie.  V'irst  40  MSEC  and  last  14  variables 
ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  this  analysis. 


RESULTS  OF  RUN:  B“3 


VARIABLES  SELECTED-  (in  m-rW) 

54,  44,  62,  38 

PERFORMANCE 

Discriminability 

Predictive 

of  Data 

Power  of  DF 

35.  90  o/o 



. o/o 



INITIAL 
fU-STAT  0.87613 


FINAL 


Approx.-  f 9.8  9643 


FINAL  CLASSIFICATION  MATRIX: 


1 

2 

3 

4 

5 

6 

7 

8 
9 
A 


40 

6 

7 

8 
6 
4 
8 
6 
4 
3 


2 1 

11  12 

12  24 

1 7 


6 

3 

2 

7 

0 

3 


6 

6 

7 

8 
2 
3 


5 3 

6 2 

7 1 

26  6 
11  13 
9 5 


9 

1 

1 

2 


2 

7 

0 

0 


3 
2 
0 

4 
1 
8 
1 
1 
1 
1 


3 

5 

3 

3 

5 

3 


6 

7 

1 


3 

0 

4 

8 

2 

10 

4 

2 

4 

3 

2 

4 

l=410nm 

6 

13 

6 

9 

4 

4 

35.9°/o  2=465nm 

3=490nm 

8 

3 

7 

4=515nm 

20 

7 

1 

5=530nm 

7 

33 

9 

6=550nm 

1 

12 

38 

7=575nm 

99 

8=595nm 

A=660nm 

9 


_BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  TESTRN 


10.  Of  TIME  POINTS-  80 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED:  2-  Central  Occlpltal/EaraiiUJKrT.s  miAl:  CP 

NUMBER  or  Cases/Names  of  Groups  Analyzed-  <’-•  indicates  group  is  being 
64/1,  64/2,  64/3,  64/4,  64/5,  tested> 

64/6,  64/7,  64/8,  64/9,  64/A 

F^EVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0,  nos 

NO.  OF  STEPS  RUN-  lo 


MEMO: 


The  first  10  variables,  i.e.  first  40  MSCE  and  last  14  variables 
t.e.  last  56  MSOE  of  the  epoch  were  not  considered  in  this  analysis. 


RESULTS  OF  RUN:  B-4 
VARIABLES  SELECTED:  (in  order) 

55,  35,  63,  38,  27,  24,  33,  43,  52,  21 


50*  10 


Predictive 
Power  of  DF 
. o/o 


INITIAL 

J-STAT  0.67121 


FINAL 

0.11323 


Approx.-  F 34. 28992 1 17.76807 


-INAE  CLASSIFICATION  MATR DC : 


0 2 4 

34  13  2 

16  33  6 

4 7 2: 

3 4 5 

6 8 II 

5 0 7 

7 4 0 

10  2 
2 0 0 


4 2 4 3 

2 3 10 

6 12  0 
23  10  4 5 

5 27  12  4 

10  15  8 6 

7 0 3 3 

0 4 3 2 

2 0 0 1 


32  8 3 2 

2 31  10  2 

1 7 41  12 

0 3 5 47 


l=410nm 
2=465nm 
, o , 3=490nm 

50-10/o  4=515nm 
5=530nm 
6=550nm 
7=575nm 
8-595nm 
9=620nm 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


DATA  CODE: 


O.  OF  TIME  POINTS:  80 


SAMPLING  RATE  4 MSEC. 


HAN'NEL  PROCESSED;  1-Top  Occipital/Ears  SUBJECT'S  INITIAL1  Cl? 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

64/1,  64/2,  64/3,  64/4,  64/5,  64/6, 

64/7,  64/8,  64/9,  64/A 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

O.  OF  STEPS  RUN:  10 


MEMO: 


RESULTS  OF  RUN:  B-5 


ARIABLES  SELEC1  (in  order) 

53,  29,  63,  32,  42,  25,  28,  23,  21,  48 


ERFORMANCE 
)iscriminability 
of  Data 

50  *78  0y/° 


Predictive 
Power  of  DF 
. o/o 


INAL  CLASSIFICATION  MATRIX- 


34  1 

7 17 


INITIAL 
-STAT  0.  70841 


FINAL 


0.19956 


Approx,-  F 28. 8121b  14.08724 


23456789A 
1 11042074 

27  16  9121230 

16  26  9120240 

10  2 35  551210 


0 50.78  /o 

3 


13  0 

16  5 6 

2 29  10 

l,n,5  49 


l=410nm  

2=465nm  Key 
3=490nm  R = Red 
4=515nm  G ~ Green 
5=530nm  B = Blue 
6=550nm  Y = Yellow 
7=575nm  w = White 
8=595nm  P = Purple 
9=620nm  BG  = Blue-Grn 
A=660nm 


BMP  07  M STEPWISE  DISCRIMINANT  ANALYSTS 


1 


* 


RUN  CODE:  TESTRN DATA  CODE: 

NO.  OF  TIME  POINTS:  80  SAMPLING  RAIL-  4 MSEC. 


CHANNEL  PROCESSED;  3-Bottom  Occipital/Ears^ju jjj FC’P'S  INITIAL1  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

64/1,  64/2,  64/3,  64/4,  64/5, 

64/6,  64/7,  64/8,  64/9,  64/A 

F-LEVEL  FOR  INCLUSION:  o.OlO  F-LEVEL  FOR  DELETION:  0. 005 

NO.  OF  STEPS  RUN:  10 


MEMO: 


The  first  10  variables,  ie.  first  40  MSEC  and  last  14  variables 
ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  this  analysis. 


Variables  selected:  (in  order) 


RESULTS  OF  RUN:B-6 


35,  54,  44,  32,  59,  36,  66,  55,  20,  46 


[performance 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

31*70  °/° 

. o/o 

- 

FINAL  CLASSIFICATION  MATRIX: 


1 

2 

3 

4 

5 

6 

7 

8 
9 
A 


46 

7 

7 

6 

4 

5 
5 

14 

5 

5 


1 2 
20  9 
10  24 

10  11 
5 10 


2 

5 

I 

3 

4 


5 

4 

4 

1 

0 


1 

2 

2 

4 

4 

5 
2 
0 
0 
0 


0 

3 

6 

13 

16 

6 

9 

6 

3 

3 


2 

5 

6 
6 
6 
13 
5 
5 
2 
3 


2 

8 

3 
2 
7 
9 

20 

4 
7 

5 


6 

4 

4 

7 

3 

3 

4 

5 
3 
5 


1 

1 

1 

1 

6 

7 

5 

16 

26 


3 

3 

3 

4 
3 
9 

5 
5 

14 


_LQ2_ 


10  29 


INITIAL 

FINAL 

U-STAT  0.87796 

0.44194 

Approx,-  F 9.72992 

6.00219 

l=410nm 

2=465nm 

Key 

3=490nm 

R = Red 

0 4=51 5 nm 

G = Green 

31.70  /o  5=530nm 

B = Blue 

6=55  Onm 

Y = Yellow 

7=575n.n 

W = White 

8=595nm 

P = Purple 

£=620nm 

BG  = Flue-Gr 

A=660nm 

_u 

BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


DATA  CODE: 


O.  OF  TIME  POINTS:  80 


SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED:  5-Right  Occipital  /Ears  SUBJECT'S  INITIAL*  CP 

"iUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

64/1,  64/2,  64/3,  64/4,  64/5,  tested) 

64/6,  64/7,  64/8,  64/9,  64/A 


IF-LEVEL  FOR  INCLUSION:  0.010 


F-LEVEL  FOR  DELETION:  0.005 


iO.  OF  STEPS  RUN:  10 


MEMO: 


The  first  10  variables,  ie.  first  40MSEC  and  last  14  variables 
ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  this  analysis. 


RESULTS  OF  RUN:  B-7 


ARIABLES  SELECTED:  (in  order) 


64,  55,  .>3,  25,  19,  44,  66,  60,  54,  27 


>ERFORMANCE 
Mscriminability 
of  Data 

34*  00  0/0 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
■STAT  0.84798 


FINAL 


0.42593 


Approx,-  F 12. 548981  6.29086 


INAL  CLASSIFICATION  MATRIX: 


_3 4 5 6 7 

I 4 0 0 2 

9 6 4 5 3 

II  11  7 6 2 

6 22  8 3 1 

10  31  5 3 

6 8 12  5 

3 13  7 1 


10  3 


13  7 
4 5 


7 8 

2 5 

3 7 
2 6 
1 2 
3 3 
5 6 
11  7 


7 10  9 

10  11  2 

5 31  10 

2 8 33 


l=410nm 

2=465nm  Key 
3=490nm  R 
4=515nm  G 


34.00  /o  5=530nm 
6=55  Onm 
7=575nm 


= Red 
= Green 
= Blue 
= Y ellow 
= White 


8=595nm  p = Purple 
9=620nm  BG  = Blue-Gi 
A=660nm 


BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


l 


RMD  07 M STEPWISE 

DISCRIMINANT  ANALYSIS 

! 

l 

RUN  CODE:  CLUSTR 

DATA  CODE: 

NO.  OF  TIME  POINTS:  80 

SAMPLING  RATE  - 4 MSEC. 

o 

CHANNEL  PROCESSED:  1 

SUBJECT'S  INITIAL:  CP 

number  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

0 

32/Red,  32/Green,  32/Blue, 

-32/1,  -32/3,  -32/5,  -32/6,  -32/  , 

F-LEVEL  FOR  INCLUSION:  0.010 

tested) 

-32/8,  -32/9 

F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  10 

0 

MEMO:  The  fir8t  10  variable8,  ie.  first  40  MSEC  and  last  14  variables 

ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  this  analysis. 


(VARIABLES  SELECTED:  (in  order) 


RESULTS  OF  RUN:  B"8 


53,  44,  62,  29,  67,  48,  38,  42,  74,  51 


(performance  1 

liscriminability 

Predictive 

of  Data 

Power  of  DF 

91*67  0//° 

rz 

o/o 

FINAL  CLASSIFICATION  MATRIX: 

IT 

R 

G B 

W~ 

31 

0 

1 

G 

0 

28 

4 

B 

0 

3 

29 

1 

1 

8 

23 

3 

1 

23 

8 

5 

0 

14 

18 

6 

2 

11 

19 

7 

3 

10 

19 

9 

15 

1 

16 

INITIAL 

FINAL  ! 

U-STAT  0.41869 

o.o8§as  I 

Approx,-  F 64.56099 

20.60997  j 

1=41  Onm 
2=465nm 
3=490nm 
4=5l5nm 
5=530nm 
6=5  5 Onm 
7=575nm 
8=595nm 
9=620nm 
A=660nm 


..1 


* 


BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  CLUSTR 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

S A M PLIN  CRATE  M SEC . 


CHANNEL  PROCESSED:  2-Central  Occipital/Ear8slJ]}jEC.|1,s  HsuriAL:  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

32/Red,  32/Green,  32/Blue,  tested) 

-32/1,  -32/3,  -32/5,  -32/6,  -32/7,  -32/8,  -32/9 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  10 

MEMO* 

" The  first  10  variables,  ie.  first  40  MSEC  and  last  14  variables 

ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  the  analysis. 


RESULTS  OF  RUN:  B-9 

VARIABLES  SELECTED:  (in  order) 

52,  73,  54,  35,  38,  63,  42,  25,  32,  26 


‘ERFORMANCE 


Diacrim  inability 
of  Data 

Predictive 
Power  of  DF 

89  * 58  0//° 

. o/o 

U 

INITIAL 

FINAL 

U-STAT  0.53898 

0.09510 

Approx.-  F 39.77399 

18.83817 

■INAL  CLASSIFICATION  MATRIX: 


1 

1 

1=41  Onm 

27 

5 

2=465nm  Key 

2 

29 

3=490nm  R 

Red 

10 

17 

4=515nm  G = 

Green 

6 

25 

5=530nm  B = 

Blue 

14 

15 

6=55  Onm  Y = 

Yellow 

14 

8 

7=57  5nm  W = 

White 

15 

12 

8=595nm  P = 

Purple 

6 

17 

9=620nm  BG  = 

Blue-G 

3 

15  105 

A=660nm 

w 


RUN  CODE;  TESTRM 


DATA  CODE- 


(NO.  OF  TIME  POINTS. 


80 


SAM  PUN  (1  HATE  4 MSEC. 


CHANNEL  PBOCESSEP:  1-Top  Oecipi.al/Ears;  SUWECT'S  INrriAl  . CP 

NUMRrn  in  / 2-Central  Occipital/Ears 
NUMBER  of  Cases/Nampa  of  Gi’ouds  Anulv/oH  n t • 

— 2L.1  louP9  Analyzed:  (’-•  indicates  group  is  being 

32/1,  32/5,  32/A  tested) 

-32/4,  -32/F,  -32/T 

F-LEVEL  FOR  INCLUSION-  0.010 


{NO.  OF  STEPS  RtlM.  10 


F-LEVEL  FOR  DELETION.  0.005 


MEMO: 


. Variables  l-io  ie.  first  40  MSEC  and 
variables  67-80  ie.  last  56MSEC  were 
not  considered  in  thisf  analysis. 

2.  Channels  1 and  2 were  simultaneously 
analyzed  by  choosing  the  odd  variables  from 
CHI  and  even  from  OH? 


[VARIABLES  SELECTED-  (in  order) 

_ Channel  2 Variables 
J ' .1  I 


38,  6g,  3^,  2^,  i,  4j),  |l,  ^ 23 

1 .t ■ i * 


Channel  1 Variables 


Performance  — I 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

93.75  0/0 

i 84*88  °/° 

ggJAL  CLASSIFICATION 


initial 

FINAL 

U-STAT  0.59704 
Approx.-  F 31. 38367 

0.08669 
20. 12962  | 

1 

5 

A 

1 

31 

0 

1 

5 

3 

29 

0 

A 

0 

2 

30 

U 

24 

3 

5 

F 

6 

26 

0 

T 

0 

1 

31 

l=Violet 
5=Green-YellowjR 
A=Red 
U= Violet 
F=Green- Y ello^JY 
T=Red  L 


G 

B 


[Key 

= Red 
= Green 
= Blue 
= Yellow 
= White 
- Purple 
fBC  = Blue-Grn 


T 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN DATA  CODE: 


NO.  OF  TIME  POINTS:  80  SAMPLING  RATE  - 4 MSEC. 

CHANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL1  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

32/1,  32/5,  32/A, 

-32/4,  -32/F,  -32/T 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0 005 

NO.  OF  STEPS  RUN:  10 


MEMO: 

Variables  1-10  ie.  first  40  MSEC  and  variables  67-80 
ie.  last  56  MSEC  were  not  considered  in  this  analysis. 


[VARIABLES  SELECTED:  (in  order) 


RESULTS  OF  RUN;  C-2 


63,  31,  52,  45,  48,  24,  21,  28,  66,  55 


[performance  I 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

. o/o 

. • o/o 

u - 

[FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.59479 

0.09828 

Approx.-  F 31. 6786( 

18.39432 

l=Violet 

5=Green- Yellow 

A=Red 

U=Violet 

F=Green-Yellow 

T=Red 


Key 

R 

G 

B 

Y 

W 

P 

IBG 


Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G: 


BMP  07  M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


O.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


HANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL*  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  bsing 

tested) 

32/1,  32/4,  32/5,  32/6,  32/A 
-32/0,  -32/C,  -32/F,  -32/S,  -32/T 

F-LEVEL  FOF  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:0* 005 

O.  OF  STEPS  RUN:  1° 


MEMO: 


Variables  1-10  ie.  first  40  MSEC  and  variables 

67-80  ie.  ^ast  56  MSEC  of  epoch  were  not  considered 
in  this  analysis. 


RESULTS  OF  RUN*.  C-3 


VARIABLES  SELECTED:  (in  order) 

63,  31,  55,  28,  42,  25,  21,  23,  37,  49 


>ERFORMANCF  

)iscriminability  Predictive 

of  Data Power  of  DF 

no  °/°  o/o 


69  38 


.66  88 


INITIAL 
I-STAT  0.62170 


FINAL 


0.12301 


Approx.-  f 23.579541  10.24572 


INAL  CLASSIFICATION  MATRIX: 


l=Violet 

4=Green 

5=Grn.  -Yellow 

6=Yellow-Grn. 

A=Red 

0=Violet 

C=Green 

F=Grn.  -Yellow  | 

S=Yellow-Grn. 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


I • 


% 


9 


9 


RUN  CODE:  TESTRN  DATA  CODE: 


NO.  OF  TIME  POINTS:  80  SAMPLING  RATE  4 MSEC. 

CHANNEL  PROCESSED:  2-Central  Occipital/Ear ECT.S  INITIAL:  CP 

NUMBER  of  Cases/Namrs  of  Groups  Analyzed:  ('-'  indicates  group  is  being 

, , tested) 

32/1,  32/4,  32/5,  32/6,  32/A 

-32/0,  -32/C,  -32/F,  -32/S,  -32/T 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  10 


MEMO: 

1.  Variables  1-10  ie.  first  40  MSEC  and  variables  67-80 

ie.  last  56  MSEC  of  the  epoch  were  not  considered  in  this  analysis. 

2.  This  run  was  made  along  with  another  similar  run  where  the 
data  used  was  preprocessed  through  a digital  filter. 


O 


o 


c 


t 


L 


RESULTS  OF  RUN: 

C-4 

VARIABLES  SELECTED:  (in  order) 

37,  54,  64, 

42,  28, 

24, 

30, 

49, 

46,  33 

PERFORMANCE 

INITIAL 

FINAL 

liscrim  inability 

I Predictive 

U-STAT  0.64535 

0.130SS 

of  Data 

1 Power  of  DF 

1 

66  * 25  0//° 

.57  -50  0/0 

Approx.-  F 21.29509 

9.85034  | 

FINAL  CLASSIFICATION  MATRIX: 

1 4 

5 

6 

i 

\ 

1 

29  1 

0 

1 

1 

l=Violet 

4 

3 19 

5 

3 

2 

4=Green 

Key 

5 

4 4 

15 

9 

) 

5=Grn. -Yellow 

R 

= Red 

6 

1 5 

7 

15 

1 

6=Yellow-Grn. 

G 

= Green 

A 

1 0 

0 

3 

28 

A=Red 

B 

= Blue 

0 

23  2 

1 

3 

3 

0=Violet 

Y 

= Yellow 

C 

1 19 

7 

2 

3 

C=Green 

W 

= White 

F 

1 8 

19 

4 

0 

F=Grn.  -Yellow 

P 

= Purple 

S 

7 8 

6 

8 

3 

S=Yellow-Grn. 

BG 

= Blue-Gr 

T 

1 0 

1 

1 

23 

T=Red 

1 1 

109 


I 


I 


l ■ 


BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


DATA  CODE: 


NO.  OF  TIME  POINTS:  80 

SAMPLING  RATE  - 4 MSEC. 

CHANNEL  PROCESSED:  2-Central  Occipital/Eara-TJDJ  ECT'S  INITIAL:  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  | 

indicates  group  is  being 

32/1,  32/4,  32/5,  32/G,  32/A, 
-32/U,  -32/C,  -32/F,  -32/S,  -32/T 

tested) 

F -LEVEL  FOR  INCLUSION:  0.010 

F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN.  10 

MEMO: 


[VARIABLES  SELECTED:  (in  order) 


[PERFORMANCE  | 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

o 

o 

o 

• 

o 

t- 

& : 0/0 

1.  Variables  1-10  ie.  first  40  MSEC  and  variables  67-80 
ie.  last  56  MSEC  were  not  considered. 

2.  The  data  was  preprocessed  through  a low-pass  digital  filter 
with  time  varying  time  constants,  which  also  removed  "DC" 
by  "zero-meaning"  it. 

3.  The  plot  of  first-2-canonical  variables  was  stratified  on  the  third  one. 

RESULTS  OF  RUN:  C-5 


55,  66,  38,  45,  28,  25,  33,  24,  22,  49 


[FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.52857 

0.07873 

Approx.-  F 34.56050 

13.25951 

A 

1 

l=Violet 

0 

A 

4=Green 

Key 

u 

5=Grn.  -Yellov 

R = Red 

3 

6= Yellow -Grn 

G = Green 

29 

A=Red 

B = Blue 

3 

0=V  iolet 

Y = Yellow 

1 

C=Green 

W = White 

0 

F=Grn-Yellow 

P = Purple 

2 

S=Yellow-Grn 

BG  = Blue-Gr 

31 

T=Red 

1 

no 


: 


BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  C1D124 
NO.  OF  TIME  POINTS:  80 


DATA  CODE: 

SAM  DUNG  RATE  4 MSEC, 


:HANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL:  CP 

DUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

30/Red,  30/Green,  30/Blue 


F-LEVEL  FOR  INCLUSION:  0.010 
NO.  OF  STEPS  RUN:  15 


F-LEVEL  FOR  DELETION:  0.005 


MEMO: 


The  results  are  reported  after  10  steps. 


RESULTS  OF  RUN:  C-6 

VARIABLES  SELECTED:  (in  order) 

74,  37,  39,  52,  18,  72,  41,  12,  68,  80 


>ERFORMANCE 

)iscriminability 

of  Data 

o/o 


Predictive 
Power  of  DF 
. o/o 


91  * 11  ' ' 1 

[FINAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

R 27 

3 

0 

G 2 

27 

1 

B 0 

2 

28 

INITIAL 

FINAL 

U-STAT  0.66125 

0.14972 

Approx,-  F 22.28453 

12.35823 

Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


[NO.  OF  TIME  POINTS:  80 


CHANNEL  PROCESSED- 1-Ton  Occipital/Ears  SUBJECT'S  INITIAL:  CP 


RUN  CODE:  GRNRAW 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


^UMBER  ol  Cases/Names  of  Groups  Analyzed:  (•-•  indicates  group  is  being 

tested) 


10/Red,  10/Green,  10/Blue,  -30/WARGRN 
[F^LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 


[NO.  OF  STEPS  RUN:  10 


MEMO: 


Red,  Green,  Blue  are  10  averages  each  computed 
from  sets  of  6 epochs  per  average. 




RESULTS  OF  RUN:  C-7 

VARIABLES  SELECTED:  (in  order) 

39,  37,  58,  75,  48,  57,  13,  12,  1,  28 


. 


PERFORMANCE  | 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

ioo*  oo  0/f° 

£3  *33  °/° 

FINAL  CLASSIFICATION  MATRIX- 


RA 

GA 

BA 

RA 

10 

0 

0 

GA 

0 

10 

0 

BA 

0 

0 

10 

WARGRN 

12 

13 

5 

INITIAL 

FINAL 

U-STAT  0.22331 

0.00076 

Approx,-  F 46.  95537 

63.  60176 

RA=Red  Average 

Key 

GA=Green  Aver. 

R = Red 

BA=Blue  Aver. 

G = Green 

WARGRN=Raw 

B = Blue 

Green 

Y = Y ellow 

W = White 

P = Purple 

BG  = Blue-G| 

112 
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BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  AVERAG 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

SAMPLING  RATE  4 MSEC. 


CHANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL;  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

10/Red,  10  Green,  10/B!ue, 

-5/Avergrn 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0. 005 

NO.  OF  STEPS  RUN:  10 

Red=  are  10  averages  of  6 epochs  each 
Green=  are  10  averages  of  6 epochs  each 
Blue=  are  10  averages  of  6 epochs  each 
AVEGRN=  are  5 averages  of  6 epochs  each 


RESULTS  OF  RUN:  C-7a 

VARIABLES  SELECTED:  (in  order) 

39,  37,  58,  75,  418,  57,  13,  12,  1,  28 


[PERFORMANCE 

Discrim  inability 
of  Data 

Predictive 
Power  of  DF 

100  *00  0^° 

0 

0 

0 

0 

• 

O 

> 

■'INAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.22331 

0.00076 

Approx.-  f 46.9553' 

63.  60176 

RED 

GREEN 

BLUE 

RED 

10 

0 

0 

GREEN 

0 

10 

0 

BLUE 

0 

0 

10 

AVEGRN 

2 

2 

1 

Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-Grn. 


BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  RAWRGB 
|no.  of  TIME  POINTS:  80 


DATA  CODE: 

SAMPLING  RATE  4 MSEC. 


:HANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT’S  INITIAL:  CP 

DUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

“ “ tested) 

30/Red,  30/Green,  30/Blue, 

-30/REDRAW,  -30/GRNRAW,  -30/BLURAW 


F-LEVEL  FOR  INCLUSION:  0.010 
NO.  OF  STEPS  RUN:  10 


F-LEVEL  FOR  DELETION:  0.C05 


MEMO: 


The  first  10  variables  ie.  the  first  40  MSEC  of  the  epoch  were 
not  considered  in  this  analysis 


RESULTS  OF  RUN:  C-8 


VARIABLES  SELECTED:  (in  order) 

74,  37,  39,  54,  72,  68,  35,  41,  19,  63 


>ER  FORMANCE 
iiscrim  inability 

of  D.  ta 


Predictive 
Power  of  DF 
_ . o/o 


QO  -00  163  * 33  — J 

FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.60871 

0.11303 

Approx,-  F 27.96211 

15.4009 

Red 
Green 
Blue 
Yellow 
White 
Purple 
Blue-Grn. 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  AVERGB 


O.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


HANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL-  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

10/REDAVE,  10/GRNAVE,  10/BLUAVE, 

-30/Red,  -30/Green,  -30/Blue 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 


O.  OF  STEPS  RUN:  10 


MEMO: 

REDAVE  = are  10  averages  of  6 epoch  each 
GRNAVE=  are  10  averages  of  6 epoch  each 
BLUAVE  = are  10  averages  of  6 epoch  each 


r 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  PA1RRB 


DATA  CODE: 


NO.  OF  TIME  POINTS-  80 


SAMPLING  RATE  -4  MSEC. 


CHANNEL  PROCESSED;  1-Top  Occipital/Ears  SUBJECT’S  INITIAL:  CP 


|NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

10/REDAVE,  10/BLUAVE,  tested) 

-10/GRNAVE,  -30/GRNRAW,  -30/REDRAW,  -30/BLURAW 

IfcM.VgL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0-000 

|hO.  OF  STEPS  RUN:  10 


MEMO: 


1.  REDAVE  = are  10  averages  of  6 epochs  each 
BLUAVE  = are  10  averages  of  6 epochs  each 
GRNAVE  = are  10  averages  of  6 epochs  each 

2.  Predictive  power  of  DF  is  computed  for  RED  and  BLUE  only. 


RESULTS  OF  RUN:  C-10 

VARIABLES  SELECTED:  (in  order) 

75,  41,  52,  35,  15,  65,  47,  48,  36,  12 


PERFORMANCE  I 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

100*00  0//° 

,96  * 67  °/° 

FINAL  CLASSIFICATION  MATRIX: 


INITIAL 

FINAL 

U-STAT  0.10192 

6.00098 

Approx,-  f 158.608] 

5 914.23120 

« 

RAVE 

BAVE 

RAVE 

10 

0 

BAVE 

0 

10 

tear 

GAVE 

8 

2 

r - 

Red 

G 

19 

11 

G = 

Green 

R 

27 

3 

B = 

Blue 

B 

5 

25 

y = 

Yellow 

w = 

White 

p = 

Purple 

116 

BG  = 

Blue-Gji 
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BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  PA1RGB 


O.  OF  TIME  POIN  TS:  80 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT’S  INITIAL: 

UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

10/GRNAVE,  10/BLUAVE, 

-10/REDAVE,  -30/REDRAW,  -30/GRNRAW,  -30/BLURAW 
F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

O.  OF  STEPS  RUN:  10 


MEMO: 

1.  GRNAVE=  are  10  averages  of  6 epochs  each. 

BLUAVE-  are  10  averages  of  6 epochs  each. 

REDAVE=  are  10  averages  of  6 epochs  each. 

2.  Predictive  power  of  DF  is  computed  for  GREEN  and  BLUE  only. 


RESULTS  OF  RUN:  C-ll 


VARIABLES  SELECTED:  (in  order) 

37,  75,  66,  51,  40,  31,  32,  24,  30,  17 


>ERFORMANCE  

)iacriminability  Predictive 

of  Data  Power  of  DF 


100*00 


85*00 


INITIAL 
-STAT  0.32701 


FINAL 
0. 00082 


Approx,-  F 37. 044911  1095.89941 


=TNAL  CLASSIFICATION  MATRIX: 


Yellow 


Blue-G 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  PA  IRE  G 
NO.  OF  TIME  POINTS;  80 


DATA  CODE; 

SAM  PUNK  RATE  4 MSEC. 


CHANNEE  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL:  CP 
number  of  Cases/Names  of  Groups  An1.lv7.rri.  indicates  group  is  being 

io/redave,  io/grnave,  ^ 

-10/BLUAVE,  -30/BLURAW , -30/REDRAW,  -30/GRNRAW 
F.LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DFT.F.TTrm.  0. Q05 

NO.  OF  STEPS  RUN:  10 


1.  REDAVE=are  10  averages  of  6 epochs  each 
GRNAVE=are  10  averages  of  6 epochs  each 
BLUAVE=are  10  averages  of  6 epochs  each 

2.  Predictive  power  of  DF  is  computed  for  RED  and  GREEN  only. 

RESULTS  OF  RUN:  C*  lU 
VARIABLES  SELECTED:  ( in  order) 

39,  14,  58,  49,  9,  53,  23,  24,  54,  26,  35 


>ER  FORM  ANC  E 
>i8criminability 
of  Data 

100*00  °/° 


Predictive 
Power  of  DF 
7R-  R7  o/o 


uo  I 76’  07  wu 

INAL  CLASSIFICATION  MATRIX- 


RAVE 

RAVE 

10 

GAVE 

0 

BLUAVE 

1 

B 

12 

R 

23 

G 

7 

GAVE 


INITIAL 
■ST AT  0.228 


FINAL 


APProx.-  F 60.  91145#  1450.26392 


Key 

R = Red 
G = Green 
B = Blue 
V = Yellow 
W = White 
P = Purple 
BG  = Blue-Grn. 


I 


BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


! t 


O 


0 


HUN  CODE:  C2SPEC 


DATA  CODE: 


NO.  OF  TIME  POINTS:  80 

SAMPLING  RATE  - 4 MSEC. 

CHANNEL  PROCESSED:  2-Central  Occipital/E arssuBJ EC T'S  INITIAL:  cp 

NUMBER  of  Cases/Names  of  Groups  Analyzed: 

indicates  group  is  being 

30/RED,  30/GREEN,  30/BLUE, 

tested) 

-30/DERRED,  -30/NRGGRN,  -30/ULBBLU 

F-LEVEL  FOR  INCLUSION:  0.010 

F-LEVEL  FOR  DELETION:0- 005 

NO.  OF  STEPS  RUN:  10 

MEMO: 


DERRED,  NRGGRN,  ULBBLU  are  simple  30  raw  epochs  of 
RED,  GREEN,  BLUE  each.  Such  1 names  were  used  to  get 
unique  letters  on  the  canonical  maps. 


RESULTS  OF  RUN:  c~12 


[VARIABLES  SELECTED:  (in  order) 


73,  57,  51,  9, 

53,  33,  47, 

58,  66,  67 

PERFORMANCE 

INITIAL 

FINAL 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

U-STAT  0.63625 
Approx,-  f 24.86957 

0.21431 

9.04907 

80  * 00  °/° 

o 

o 

o 

o 

FINAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

R 

23 

5 

2 

Key 

G 

7 

22 

1 

R = Red 

1 

B 

2 

1 

27 

G = Green 

R 

14 

6 

10 

B = Blue 

G 

6 

16 

8 

Y = Yellow 

B 

4 

11 

15 

119 

W = White 
P = Purple 
BG  = Blue-Gi 

r* 

BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


DATA  CODE: 


O,  OF  TIME  POINTS:  80 


SAMPLING  RATE  --  4 MSEC, 


HANNEL  PROCESSED:  2-Central  Occipital/Ear^mjj EC'T'S  INITIAL1 

1 UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

64/Yellow,  64/  DUMMY 


F-LEVEL  FOR  INCLUSION:  0.010 


F-LEVEL  FOR  DELETION:®- 005 


O.  OF  STEPS  RUN:  10 


MEMO: 


1.  The  first  10  variables  ie.  first  40  MSEC  and  last  14  variables 
ie.  last  56  MSEC  of  epoch  were  not  considered  in  this  analysis. 

2.  "DUMMY"  consists  of  64  epochs  of  a constant  value  of  1.0  each  at 
time  point. 


RESULTS  OF  RUN:  C-13 


ARIABLES  SELECTED:  (in  order) 


25,  62,  16,  43,  60,  48,  34,  41,  37,  44 


PERFORMANCE 
li  sc  rim  inability 
of  Data 

98  *44  o/° 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
-STAT  0.38685 


FINAL 


Approx.-  f 190. 711619  145.13947 


INAL  CLASSIFICATION  MATRIX: 


pUMMY 


Y ellow 


Blue-G 


BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


O.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


HANNEL  PROCESSED:  2-Central  Occipital/Ear kcT'S  INITIAL:  CP 


UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

64/RED,  64/DUMMY 


F-LEVEL  FOR  INCLUSION:  0.010 


O.  OF  STEPS  RUN:  10 


F-LEVEL  FOR  DELETION:  0.005 


MEMO: 


1.  The  first  10  variables  ie.  first  40  MSEC  and  the  last  14  variables 
ie.  last  56  MSEC  were  not  considered  in  this  analysis. 

2.  ’DUMMY*  consists  of  64  epochs  of  a constant  value  of  1.0  eaeh  at 
each  time  point. 


RESULTS  OF  RUN:  C-13a 


VARIABLES  SELECTED:  (in  order) 

14,  20,  22,  28,  44,  25,  55,  63,  18,  15 


>ERFORMANCT 
)iscriminability 
of  Data 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
-STAT  0.28878 


i FINAL 


0.05820 


100  00  L 

If INAL  CLASSIFICATION  MATRIX: 


Approx,-  F 310. 3198*  189.34644 


DUMMY 


R 

64 

0 

DUMMY 

0 

64 

121 

Yellow 

Whit 

e 

Purp 

>le 

Blue 

-Gi 

■ 

E 


i 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


DATA  CODE: 


NO.  OF  TIME  POINTS:  80 


SAMPLING  RATE  4 MSEC, 


CHANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

NUMBER  of  Caset /Names  of  Groups  Analyzed:  ('-•  indicates  group  is  being 

tested) 

64/GRN,  64/Dummy 


f-level  FOR  INCLUSION.  0.010 


F-LEVEL  FOR  DELETION:  n 005 


M>.  OF  STEPS  RUN- 


MEMO: 


1.  The  first  10  variables,  i.e.  first  40  msec  and  the  last 
14  variables,  i.e.  last  56  msec  of  the  epoch  were  not 
considered  in  this  analysis.  


2.  'Dummy'  consists  of  64  epochs  of  a constant  value  of 
1.0  at  each  time  point 


RESULTS  OF  RUN:  C-14 


VARIABLES  SELECTED;  ( in  order) 

14,  17,  38,  34,  43,  50,  29,  31,  16,  64 


PERFORMANCE 
Rscrimlnability 
of  Data 

100  * 00  °/° 


Predictive 
Power  of  DF 
. o/o 


INITIAL 
U-STAT 
0.42088 
Approx.-  f 


FINAL 


0.09659 


FINAL  CLASSIFICATION  MATRIX; 


py 

R = Red 
G = Green 
B = Blue 
Y = Yellow 
W = white 
P = Purple 
BG  = Blue-Grn. 


P 


KAl 


BMD  07  M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 


O.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC 


H ANN  EL  PROCESSED-.  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

1 UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 


64/BLUE,  64/Dummy 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  10 


MEMO:  l.  The  first  10  variables,  i.e.  first  40  msec  and  the  last  14  variables 
i.e.  last  56  msec  were  not  considered  in  this  analysis. 


2.  'Dummy'  consists  of  64  epochs  of  a constant  value  of  1.0  at  each  time 
point. 


RESULTS  OF  RUN:  C-16 


ARIABLES  SELECTED:  (in  order) 


53,  16,  63,  21  , 

13,  26,  27,  58, 

PERFORMANCE 

Discriminability 

Predictive 

of  Data 

Power  of  DF 

99  # 22 

- -4- 

o/o 

INITIAL 

FINAL 

U-STAT 

0.39659 

WW 

0.07121 

152.59428 

FINAL  CLASSIFICATION  MATRIX: 

_B Dummy 

B 63  1 

Dummy  0 64 


Key 
R = 

r 

Red 

G = 

Green 

B = 

Blue 

Y = 

Y ellow 

W = 

White 

P = 

Purple 

BG  = 

Blue-G^ 

l 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  TESTRN 

NO.  OF  TIME  POINTS: 


DATA  CODE; 

SAMPLING  RATE  4 MSEC. 


:HANNEL  PROCESSED:  2-Central  SU13JKC  T'S  INITIAL:  CP 

Occipital/Ears 

UJMBER  of  Cases/Naaies  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

64/ORANGE,  64/ Dummy 


F-LEVEL  FOR  INCLUSION:  0.010 
NO.  OF  STEPS  RUN;  10 


F-LEVEL  FOR  DELETION:  0.005 


— EM-°:  1 • The  first  variables,  i.e.  first  40  msec  and  the  last  14  variables, 
i.e.  last  56  msec  were  not  considered  in  this  analysis. 

2.  'Dummy1  consists  of  64  epochs  of  a constant  value  of  1.0  at  each 
time  point. 


RESULTS  OF  RUN:  C-18 

VARIABLES  SELECTED:  (in  order) 

35,  59,  12,  42,  64,  32,  13,  61  , 20,  14 


PERFORMANCE 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

98  • 44  °/° 

. o/o 

A ! 

=TNAL  CLASSIFICATION  MATRIX; 


0 

Dummy 


INITIAL 

U-STAT 

0.43029 

^ra?76F 


FINAL 

0.11126 

93.45956 


Key 

R = Red 
G = Green 
B = Blue 
Y = Yellow 
W = White 
P = Purple 
BG  = Blue-Grn. 
3 = Oranqe 


w 


BMP  07M  STEPWISE  DISCRIMINANT  ANAL Y SIS 


0 


o 


0 


RUN  CODE:  TESTRN 


NO.  OF  TIME  POINTS- 

un  in  vturr, 

80  SAMPLING  RATE  4 

MSEC. 

CHANNEL  PROCESSED- 

^tlal/Ears  S,,,iJEC  I"s  INITtA^ 

CP 

NUMBER  of  Cases/Names 

of  Gioups  Aralyzed:  1 indicates  proun  is 

tested) 

64/VIOLET,  64/Dummy 

F-LEVEL  FOR  INCLUSION- 

°-010  F-LEVEL  FOR  DELETION.  0.005 

NO.  OF  STEPS  RUN: 

10 

MEMO: 


' 1. 


Ihe  f]rsJ;  var1ables*  i-e-  first  40  msec  and  the  last  14  variables 
i.e.  last  56  msec  were  not  considered  in  this  analysis. 

'Dummy'  consists  of  64  epochs  of  a constant  value  of  1.0  at  each 
time  point. 


O 


0 


[VARIABLES  SELECTED:  (in  order) 


RESULTS  OF  RUN:  C-19 


13,  35,  18,  46,  32,  11,  73,  14,  48,  65 


PERFORMANCE 

Discrim  inability 
of  Data 

Predictive 
Power  of  DF 

96  * 88  °/o 

. o/o 

-4- 

(FINAL  CLASSIFICATION  MATRIX: 


Dummy 


V 

| Dummy 


60 

0 


4 

64 


INITIAL 

FINAL 

U-STAT 

0.46307 

0.13690 

Approx,-  f 
1—146,00867 

71  7C1Q/I 

Key 

IR  = Red 

G = Green 
B = Blue 
Y = Yellow 
W = White 
P = Purple 
BG  = Blue-Gyn. 

-V 1 0 1 < 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  TESTRN 

NO.  OF  TIME  POINTS:  80 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


CHANNEL  PROCESSED:  .SUHJEC  rS  miTIAL:  CP 

NUMBER  of  Cases/Names  of  Group.  Araftacd:  indicates  group  is  being 

tested) 

64/ VIOLET,  64/Dummy 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 


NO.  OF  STEPS  RUN: 


Blue-Grn. 


BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  PGGRBL 


O.  OF  TIME  POINTS: 


DATA  CODE:  * 


SAMPLING  RATE  - 4 MSEC. 


HANNEL  PROCESSED:  i -Top  Occipital/Ea^UJ ECT'S  INITIAL:  CP 


1UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

32/RED,  32/GREEN,  32/BLUE 


F-LEVEL  FOR  INCLUSION:  0.010 


O.  OF  STEPS  RUN:  10 


F-LEVEL  FOR  DELETION:  °-005 


MEMO:  This  is  the  oth  iteration  of  a series  of  3 iterations  where  each 
iteration  was  run  after  rejectinq  outliers. 


Outliers  to  be  rejected:  R9 

G14,  G18,  G21 , G29 
B6,  B14,  B27 


RESULTS  OF  RUN:  D-l 

VARIABLES  SELECTED:  (in  order) 


53,  44,  62,  29,  67,  48,  38,  42,  74,  51 


PERFORMANCE 

Discriminability 

Predictive 

of  Data 

Power  of  D 

91 .67  o/o 

. o/o 

A 

FINAL  CLASSIFICATION  MATRIX: 

INITIAL 

FINAL 

U-STAT 

0.41869 

Approx.-  F 
64,  56099 

0.08384 
. 20 .,,60997 

Key 

R - Red 
G = Green 
B = Blue 
Y = Yellow 
W = White 
P = Purple 
BG  ~ Blue— G 


f 


i 


DMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  STRIPD DATA  CODE; 

NO.  OF  TIM E POINTS:  80  SAMPLING  RATE  - 4 MSEC. 

CHANNEL,  PROCESSED:  1-Top  Occ i pi ta  1 /Ea |3j KO'U'S  INITIAL*  CP 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

31/RED,  28/GREEN,  29/BLUE 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  10 

MEMO:  This  is  the  first  iteration  after  rejectinq  outliers  identified  by  the 
0th  iteration. 

Outliers  to  be  rejected:  R'8,  R ' 29 

G * 1 7 , G' 28 


RESULTS  OF  RUN:  D_2 


Variables  selected:  (in  order) 


53,  44,  74,  48,  67,  62,  38,  42,  51,  33 


PERFORMANCE  | 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

95.46  °/° 

. o/o 

-4 

FINAL  CLASSIFICATION  MATRIX; 

RGB 
R 29  0 2 

G 0 26  2 

B 0 0 29 


INITIAL 

FINAL 

U-STAT 

0.37660 

0.03784 

31.47003 

Key 

R = Red 
G = Green 
B = Blue 
Y = Y ellow 
W = White 
P = Purple 
BG  = Blue-Gi 


) 


:>  j 

j 


I 


> 


) 


f 


128 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE: 


JO.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROP KSSED-  1 -TopOccipi tal/Ear*„uj ECT'S  INITIAL:  CP 

NUMBER  ol  Cases/Names  of  Groups  Analyzed:  Indicates  group  is  being 

tested) 

29/RED,  26/GREEN,  29/BLUE 

F^LEVEL  FOR  INCLUSION-  A F-T.fvft  iron  nr,  


0.010 


F-LEVEL  FOR  DELETION:  Q , 


JO.  OF  STEPS  RUN; 


MEMO;  This  is  the  second  iteration  after  rejectlnq  outliers  identified  by  the 
first  iteration.  J 


RESULTS  OF  RUN:  [)-3 


VARIABLES  SELECTED:  (In  order) 

53,  44,  74,  48,  42,  38,  67,  51,  62,  33 


*ER  FORMANC  E 
Mscriminability 

of  Data 

1 00*  00  % 


Predictive 
Power  of  DF 
. o/o 


INITIAL 

U-STAT 

0.35824 

APP^552^9 


FINAL 

0.02612 


37.35187 


•INAL  CLASSIFICATION  MATRIX: 


0 26 
0 0 


Key 

R = Red 
G = Green 
B = Blue 
Y = Y ellow 
W = White 
P = Purple 
BG  = Blue-Grn. 


BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  HAFRGB 


O.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  4 MSEC. 


IHANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

JUMPER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 


32/RED,  32/GREEN,  32/BLUE,  32/YELLOW , 

-50/1  RED,  -50/2GREEN,  -50/3BLUE,  -50/4YELL0W 


r -LEVEL  FOR  INCLUSION: 
NO.  OF  STEPS  RUN: 


0.010 

6 


F-LEVEL  FOR  DELETION:  Q 0Q5 


MEMO:  RED>  GREENj  BLUE,  YELLOW  are  from  10  Color  Experiment 

1RED,  2GREEN,  3BLUE , 4YELL0W  are  from  Opponent  Color  Experiment 


RESULTS  OF  RUN:  D-4 


MlRIABLES  SELECTED:  (in  order) 


52,  20,  76,  28,  63,  25 


PERFORMANCE  1 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

70  • 31°/° 

40  • 50  °/° 

4 — 

INITIAL 

FINAL 

U-STAT 

0.61527 

0.16216 

ApP^835f9 

16.90039 

■INAL  CLASSIFICATION  MATRIX: 
R G B 


30 

2 

0 

4 

12 

8 

1 

7 

22 

0 

3 

3 

21 

11 

9 

11 

13 

22 

10 

22 

16 

3 

9 

7 

Key 

R = Red 
G = Green 
B = Blue 
Y = Y ellow 
W = White 
P = Purple 
BG  = Blue-G 


k 


i 


BMD07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE: 

(NO.  OF  TIME  POINTS*  80 


DATA  CODE; 

SAMPLING  RATE  - 4 msec. 


PHANNEL  PROCESSED: 


^~0cc1pftal/Ears  GROT'S  INITIAL:  CP 


p UMBER  of  Caaes/Namns  nf  n,.„upe  An,lvznd.  indlcates 

tested) 

50/RED,  50/GREEN,  50/BLUE 
F-LEVEL  FOR  INCLUSION-  q. 010 

(NO.  OF  STE PS  R UN-  6 


group  is  being 


F-LEVEL  FOR  DELETION-  0.005 


resui; 

VARIABLES  SEI  ECTED-  (in  order) 


53,  18, 

25,  36,  33,  22 

PERFORMANCE  I 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

86  • 00  °/° 

• o/o 

FINAL  CLASSIFICATION  matpty. 

R G 

B 

R 45 

2 3 

G 6 40  4 

B 1 

5 44 

OF  RUN:  D-5 


INITIAL 

FINAL 

U-STAT 

0.59632 

0.21241 

27.68384 

Key 

R 

Red 

G = 

Green 

B = 

Blue 

Y = 

Yellow 

W = 

White 

P = 

Purple 

BG  = 

Blue-G: 

BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  pmrwHI 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

SAMPLING  RATE  - ^ MSEC. 


CHANNEL  PROCESSED:  2-Central  SUBJECT’S  INITIAL:  CP 

Occipital/Ears 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/PUR001  , 25/ WH 1 001, 

-25/PUR002.-25/WHI002 

F-LEVEL.  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  6 


MEMO: 


RESULTS  OF  RUN: 


ARIABLES  SELECTED:  (in  order) 

25,  75,  55,  24,  63,  72 


Discriminability 
of  Data 

Predictive 
Power  of  DF 

98  • 00  °/° 

86  • 00  °/° 

A 

INAL  CLASSIFICATION  MATRIX: 
PI  W1 


INITIAL 

U-STAT 

0.44691 


FINAL 

0.16066 

37.44101 
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BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  YELWHI 


O.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


HANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL: 

Occipital/Ears 

DUMBER  of  Cases/Names  of  Groups  Analyzed:  ( *— 1 indicates  group  is  being 

tested) 

25/YEL001 , 25/WHI001 , 

-25/YEL002 , -25/WHI002  • . nnc 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION^  0.005 

NO.  OF  STEPS  RUN:  6 


MEMO: 


RESULTS  OF  RUN: 


ARIABLES  SELECTED:  (in  orderi 


63,  65,  72,  70,  56,  38 


PERFORMANCE  

Discrim  inability  Predictive 

of  Data Power  of  DF 

94  00P/o  74  • 00  °/° 

- . 

FINAL  CLASSIFICATION  MATRIX: 


Y1  22 


Y2  14 


INITIAL 

J-STAT 

0.37577 


FINAL 


0.21019 


Appr>:7f79o L 26.92953 


Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G 


- — 


■ 


BMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  BLUWHI 


O.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  - 4 MSEC. 


HANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/BLU001 , 25/WHiOOl 
-25/ BLU002 , -25/WdI002 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.00' 

NO.  OF  STEPS  RUN:  6 


MEMO: 


RESULTS  OF  RUN:  D-8 


VARIABLES  SELECTED:  (in  order) 
54,  64,  22,  61,  33,  75 


INITIAL 

U-STAT 

0.39183 


PERFORMANCE 
Discriminability 
of  Data 

Predictive 
Power  of  DF 

100  • 00  °/° 

92  • 00  °/° 

A ■ . ■■ 

FINAL 

0.14790 

41.28970 


■INAL  CLASSIFICATION  MATRIX: 


LI 


Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G 


. ^ - 


- - 


BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  GRNWHI 


O.  OF  TIME  POINTS: 


DATA  CODE: 


SAMPLING  RATE  - A MSEC. 


HANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital /Ears 

UMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/GRN001 , 25/WHI001, 

-25/GRN002,  -25/WHI002 

-LEVEL  FOR  INCLUSION:  F-LEVEL  FOR  DELETION:  q q05 


F-LEVEL  FOR  INCLUSION: 
NO.  OF  STEPS  RUN; 


MEMO: 


0.010 

6 


RESULTS  OF  RUN;  D-9 


ARIABLES  SELECTED:  (in  order) 

66,  36,  57,  26,  68,  64 


INITIAL 

J-STAT 

0.46778 


PERFORMANCE 

Discriminability 

Predictive 

of  Data 

I-  ower  of  DF 

100*00  o/o 

k 30 

cr> 

• 

o 

o 

o 

o' 

'INAL  CLASSIFICATION  MATRIX: 

G1  W1 


Appax6i24 


FINAL 

0.17890 

32.89241 


Key 

R = Red 
G = Green 
B = Blue 
Y - Yellow 
W = White 
P = Purple 
BG  = Blue-G 
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BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CCDE;  REDWHI  DATA  CODE: 

NO.  OF  TIME  POINTS:  80  SAMPLING  RATE  - 4 MSEC. 

CHANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/RED001 , 25/WHI001 , 

-25/RED002,  -25/WHI002 

F-LEVEL  FOR  INCLUSION:  Q F-LEVEL  FOR  DELETION:  Q 0Q5 

NO.  OF  STEPS  FUN:  #• 


MEMO: 


RESULTS  OF  RUN; 


VARIABLES  SELECTED:  (in  order) 
32,  9,  42,  49,  72,  68 


PERFORMANCE 

INITIAL 

FINAL 

Di8crim  Inability 

Predictive 

U-STAT 

of  Data 

Power  of  DF 

0.39981 

0.14562 

100-  00  o/o 

90  • 00  °/° 

Appro^:nfMa 

42.04842 

FINAL  CLASSIFICATION  MATRIX: 

R1 

W1 

R1  25 

0 

Key 

W1  0 

25 

R = Red 

G = Green 

R2  21 

4 

B = Blue 

Y = Yellow 

W2  1 

24 

W = White 

P = Purple 

136 

BG  = Blue-Gi 
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BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  REDBGR 
NO.  OF  TIME  POINTS: 


DATA  CODE: 

SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  ( *— * indicates  group  is  bein^ 

tested) 

25/RED001 , 25/BGR001 

-25/RED002,  -25/BGR002 

F-LEVEL  FOR  INCLUSION:  q q-jq  F-LEVEL  FOR  DELETION:  q.005 

NO.  OF  STEPS  RUN:  6 


MEMO: 


RESULTS  OF  RUN:  D-ll 

VARIABLES  SEI.ECTED:  (in  order) 

52,  38,  23,  79,  56,  36 


>ERFORMANCE  

Hscrim  inability  Predictive 

of  Data Power  of  DF 

100.00  0/0  88  .00  0/0 


INITIAL 

FINAL 

U-STAT 

0.18705 

0.58376 

ATO2!T3f 

31.14841 

•INAL  CLASSIFICATION  MATRIX: 


Red 

Green 

Blue 

Y ellow 

White 

Purple 

Blue-G 


BMP  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  BLUYFL 
NO.  OF  TIME  POINTS: 


DATA  CODE; 

SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL-  CP 

Occipital/Ears 

NUMBER  of  Caaes/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/BLU001 , 25/YEL001 , 

-25/BLU002,  -25/YEL002 

F-LEVEL  FOR  INCLUSION:  F-LEVEL  FOR  DELETION: 

0-010  ” 0.005 

NO.  OF  STEPS  RUN:  g 


MEMO: 


RESULTS  OF  RUN:  D-12 

VARIABLES  SELECTED:  (in  order) 


33,  25,  38,  79,  24,  80 


[PERFORMANCE 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

98  . 00 0/0 

86  . 00  0/0 

-4 

INITIAL 

J-STAT 


0.46641 

ApprSS:<n<r7i 


=TNAL  CLASSIFICATION  MATRIX: 


FINAL 


0.20416 

27.93663 


Key 

R = Red 
G = Green 
B = Blue 
Y = Y ellow 
W = White 
P = Purple 
BG  = Blue-C 


I 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  MTPinu 


DATA  CODE: 


[NO.  OF  TIME  POINTS: 


80 


SAMPLING  RATE  - 4 MSEC. 


PHANNEL  PROCESSED:  2-Central  SUBJECT’S  INITIAL:  CP 

Occipital/Ears 

[NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 


50/RED,  50/GREEN,  50/BLUE, 
-50/RED003,  -50/GRN003,  -50/BLU003, 
-50^/ RED  050,  -50/GRN050,  -50/BLU050 


E-LEVEL  FOR  INCLUSION: 


tested) 


F-LEVEL  FOR  DELETION: 


[NO.  OF  STEPS  RUN: 


0.010 

6 


MEMO: 


0.003 


A Discriminant  Function  was  trained  on  higher  intensity  level  data 
(attenuation  0.0)  and  applied  on  lower  intensity  level  data 
(attenuation  0.30)  and  lower  intensity  level  data  (attenuation  0.50) 


© 

O 


RESULTS  OF  RUN:  D"13 

VARIABLES  SELECTED:  (in  order) 

53,  18, 

25,  36,  33,  22 

[PERFORMANCE 

INITIAL 

FINAL 

Discriminabilitv 

Predictive 

U-STAT 

of  Data 

Power  of  DF 

0.59632 

0.21241 

86  - 00°/o 

42  .-33  °/° 

Apr«S67f 

27.68384  1 

[FINAL  CLASSIFICATION  MATRIX: 


R 

G 

R 

45 

6 

G 

2 

40 

B 

3 

4 

86% 

Key 

R = Red 

B 

REDO. 3 
GRN0.3 

1 

38 

12 

5 

— ro — 

19 

44 

2 

19 

55% 

G = Green 

B = Blue 

BLU0.3 

9 

16 

25 

Y = Yellow 

REDO. 5 

~TJ — 

23 

“IT 

W = White 

GRN0.5 

24 

9 

17 

30% 

1 on 

P = Purple 

BLU0.5 

15 

26 

0 

BG  = Blue-G] 

: 
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BMS  °™  STEPWISE  DISCRIMTMAht  slg 


RUN  CODE:  HIGL03 

NO.  OF  TIME  POINTS:  80 


DATA  CODE; 

SAMPLING  RATE  4 MSEC. 


CHANNEL  PROCESSED;  2-Centra  1 

~ - ^ Occipital/Ears  «m»L:  CP 

— °f  Cases/Names  of  Group,  An,ly,,rt:  ln(Jicalc5  group  „ 

50/ RED003,  50/GRN003,  50/BLU003  tested) 

-50/RED000,  -50/GRN000,  -50/BLU000,’ 
r-hMM  ag-  -50/BLU050 

— ■ o 010  F-LEVEL  FOR  DELETION. 

NO.  OF  STEPS  RUN-  * “0.005 

6 


RESULTS  OF  RUN- 


VARIABLES  SELECTFn.  { in  order) 
27,  53,  66,  33,  70,  24 


Performance  — 

jTscrim  inability  Predictive 

°1  Data Power  of  DF 

83  • 00  °/°  47  . on  o/o 


INITIAL 

U-STAT 

0.69650 


FINAL 

0.24789 

3,86777 


j^AL  CLASSIFICATION  MATnrv. 


pir 

R = 

Red 

G = 

Green 

B = 

Blue 

Y - 

Yellow 

W = 

White 

P = 

Purple 

BG  = 

Blue-Gjn. 

1 


BMP  07M  STEPWISE  DISCRIMm.mT  aNALYS1S 


RUNCODE:  HIGL05 


DATA  COOK: 


<0.  OF  TIME  POTTSJ TS: 


SAMPLING  RATE  - 4 MSEC. 


CHANNEL  PROCESSED. 

2-Central 


- » ••'-'V  UUUfjU  - ^ .... 

- 2-Central  SUBJECT'S  INITIAL-  cp 

^ ~ «.-■  indicates  groUp  „ bel„B 

50/RED050,  50/GRN050,  50/BLU050  tested) 

-50/RED000,  -50/GRN000,  -SO/BLUOOO; 

p^agggr  • -w^uoao 

“ 0 010  ZrLEVEL  FOR  DELETION; 

pjQ.  OF  STEPS  RUN.  * 0.005 


MEMO:  . n.  . . I 

intensity  lK.ld.S  o?  aul^M^n  0 3M  ° a“a™a«o"  0.000  and  on 


RESULTS  OF  RUN:  0-15 

VARIABLES  SELECTED  (in  order) 


30,  66,  59,  22,  33,  36 


»ER  FORMA  NCE 
Hscriminability 
of  Data 


56  • 67 0/0 


Predictive 
Power  of  DF 

i • ««  % 


INITIAL 
J-STAT  ~ 


FINAL 


0.70636 


0.24912 

23.75012 


jNAL  CLASSIFICATION  MATRIX- 


R0.5 
GO.  5 
B0. 5 
RO.O 
GO.  0 
B0. 0 
R0.3 
GO. 3 
B0. 3 


Key 

R = Red 
G = Green 
B = Blue 
Y = Yellow 
W = White 
P = Purple 
BG  = Blue-Grn. 


r 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE: DATA  CODE: 

NOOF  TIME  POINTS:  80  SAMPLING  RATE  , 4 MSFr 

— -L  PHOCKSSER  2_0cc1pUd]/Ears  mmAL:  CF 

^UMBER  of  Cases/Names  of  Groups  Analyzed:  (•-•  indicates  group  is  being 

50/RED,  50/GREEN,  50/BLUE  tG8ted) 

-100/REDN0N,  -100/GRNN0N,  -100/BLUNON 

F-LEVEL  FOR  INCLUSION-  A A1A  p.t.vvpt  


0.010 


F-LEVEL  FOR  DELETION* 


O.  OF  STEPS  RUN: 


= 50  PURPLE  + 50  YELLOW 
* 50  BLU-GRN  + 50  YELLOW 
= 50  BLU-GRN  + 50  PURPLE 


RESULTS  OF  RUN;  D-17 
VARIABLES  SEI  ECTED:  (in  order) 

53,  13,  25,  36,  33,  22 


PERFORMANCE  

Mscriminability  Predictive 

of  Data Power  of  DF 

86  • 00%  33  - 67  o/o 


jfflAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

REDN0N 

GRNN0N 

BLUN0N 


R 

G 

B 

45 

2 

3 

6 

40 

4 

1 

5 

44 

28 

25 

47 

24 

10 

66 

12 

25 

63 

INITIAL 

FINAL 

U-STAT 

0.59632 

ApP^756f4 

0.21241 

27.68384 

Key 

R = Red 
G = Green 
B = Blue 
Y = Yellow 
W = white 
P = Purple 
BG  = B1 


BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  APYHAF 


O.  OF  TIME  POINTS:  80 


DATA  CODE:  


SAMPLING  RATE  - 4MSEC, 


HANNEL  PROCESSED:  2-Central  SUBJECT'S  INITIAL:  CP 

Occipital/Ears 

NUMBER  of  Cases/Names  of  Groups  Analyzed:  indicates  group  is  being 

tested) 

25/BGR,  25/  PURPLE,  25/YELLOW 

-25/BGR002,  -25/PUR002,  -25/YEL002 

F-LEVEL  FOR  INCLUSION:  0.010  F-LEVEL  FOR  DELETION:  0.005 

NO.  OF  STEPS  RUN:  6 

MEMO: 


RESULTS  OF  RUN:  Q-l 9 

VARIABLES  SELECTED:  (in  order) 

25,  33,  66,  47,  26,  29 


PERFORMANCE  | 

Mscrim  Inability 

Predictive 

of  Data 

Power  of  DF 

90  • 67  °/° 

73  • 33  % 

■INAL  CLASSIFICATION  MATRIX: 


BG 

P 

Y 

BG 

22 

2 

1 

P 

2 

23 

0 

Y 

2 

0 

23 

BG2 

17 

3 

5 

P2 

10 

13 

2 

Y2 

0 

0 

25 

INITIAL 

FINAL 

U-STAT 

0.42022 

ApTO654f 

0.09934 

24.26294 

Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G 
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BMD  07 M STEPWISE  DISCRIMINANT  ANALYSIS 


'1 


RUN  CODE:  rrrpri DATA  CODE: 

NO.  OF  TIME  POINTS:  8Q  SAMPLING  RATE  = ^ MSEC. 

CHANNEL  PROCESSED:  1-Top  Occ 'Dital/Ears  SUBJECT'S  INITIAL:  CP 

NUMBER/NAMES  OF  GROUPS  ANALYZED: 

32/RED,  32/GREEN,  32/BLUE 

F-LEVEL  FOR  INCLUSION;  0.010  F-LEVEL  FOR  DELE  TION:  0.QQ5 

NO.  OF  STEPS  RUN:  10 

MEMO: 


ARIABLES  SELECTED:  (in  order) 


53,  44,  62,  29,  67,  4B,  38,  42,  74,  51 


PERFORMANCE 

Discriminability 

of  Data 

91..  67  o/o 


Predictive 
Power  of  DF 
. o/o 


INITIAL 

U-STAT 

0.41869 

AppWMfrF 


FINAL 


0.08384 

20.60997 


FINAL  CLASSIFICATION  MATRIX: 


Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G] 


BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE;  RRRRR! DATA  CODE; 

NO.  OF  TIME  POINTS;  80  SAMPLING  RATE  = 4'  MSEC. 


CHANNEL  PROCESSED:  1-Top  Occipital/Ears  SUBJECT'S  INITIAL:  CP 

NUMBER/NAMFS  OF  GROUPS  ANALYZED: 

32/RED,  32/GREEN,  32/BLUE, 

-32/A,  -32/4,  -32/2 

F-LEVEL  FOR  INCLUSION;  Q Q1g  F-LEVEL  FOR  DELE  TION:  Q 0Q5 

NO.  OF  STEPS  RUN:  10 


MEMO:  jhe  first  io  variables,  i.e.  the  first  40  msec  of  the  epoch  are  nfft 
$ considered  in  this  analysis 


0 


c 


c 


RESULTS  OF  RUN:  D_2i 


Variables  selected:  (in  order) 


75,  63,  53,  42,  21,  59,  55,  38,  35,  25 


PERFORMANCE 

Discriminability 
of  Data 

Predictive 
Power  of  DF 

91  * 67°/° 

^4 

• 

o 

00 

o 

> 

[FINAL  CLASSIFICATION  MATRIX: 


R 

G 

B 

R 

31 

0 

1 

G 

1 

27 

4 

B 

0 

2 

30 

A 

26 

1 

5 

4 

0 

28 

4 

2 

1 

11 

20 

145 


INITIAL 


IU-STAT 

0.58915 
roxin,ate-F 
P-423R5 


FINAL 


0.12531 

L32954- 


Key 

R 

G a 
B = 
Y = 
W = 
P = 
IBG  = 


Red 
Green 
Blue 
Yellow 
White 
Purple 
Blue-Grn. 


r 


HMD  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


I 


RUN  CODE: 

*NO.  OF  TIME  POINTS: 


DATA  CODE: 


80 


SAMPLING  RATE  = , 4 MSEC. 


[channel  PROCESSED:  2-Central  SUBJECT'S  INITIAI  • TP 

, Occipital/Ears  ’ 

PLUMBER /NAMES  OF  GROUPS  ANALYZED; 

25/RED,  -25/1  RED,  25/GRN,  -25/2GRN,  25/BLU , -25/3BLU 


FiLEJ/FL  FOR  INCLUSION:  0.010 
[NO.  OF  STEPS  RUN-  6 


F-LEVEL  FOR  DELE  TION:  O.OoJ 


MEMO: 


OnW  the  foil  owl  r,  i windows  were  considered  in  this  run- 
JJT  arla  es  2-®TTi7  92  msec-116msec, 

2 ara  es  34-38,  i.e.  136  msec-152  msec, 

a uar  3k  65  i-e>  172  msec-l 92  msec, 

(4)  Variables  51-56,  I.e.  204  msec-234  msec. 


RESULTS  OF  RUN:  D-23 


[VARIABLES  SELECTED:  (in  order) 

53,  23,  36,  48,  43,  25 


PERFORMANCE  I 

DIscriminability 
of  Data 

Predictive 
Power  of  DF 

78*67  °/° 

■“CTl 

• 

o 

o 

o 

> 

INITIAL 

FINAL 

U-STAT 

0. 56367 
Approximate-F 
21  86734 

0.26751 
■iq  .42.31 7 

INAL  CLASSIFICATION  MATRIX- 


R 

G 

B 

1R 

2G 

3B 


21 


1 

19 

5 

Key 

1 

5 

19 

R 

Red 

17 

•> 

C 

6 

G = 

Green 

2 

13 

10 

B = 

Blue 

0 

7 

18 

Y = 

Yellow 

W = 

White 

P = 

Purple 
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BG  = 

Blue-G: 

BMP  07M  STEPWISE  DISCRIMINANT  ANALYSIS 


RUN  CODE:  GRNP1IR  

NO.  OF  TIME  POINTS:  80 

CHANNEL  PROCESSED:  2.Central 

, Occipital /Ears 

NUMBER /NAMES  OF  GROUPS  ANALYZED: 


DATA  CODE: 

SAMPLING  RATE  = 4 MSEC. 

SUBJECT'S  INITIAL:  CP 


25/GRN001,  25/PUR001. 
25/GRN002,  -25/PUR002 


F-LEVKL  FOR  INCLUSION:  Q 01Q 
NO.  OF  STEPS  RUN;  6 

MEMO: 


RESULTS  OF  RUN:  q_22 


VARIABLES  SELECTED;  (in  order) 
66,  36,  15,  25,  29,  19 


PERFORMANCE 
Discriminability  Predictive 
of  Data Power  of  DF 

98  • 00°/°  | 78-  00  °/° 

"INAL  CLASSIFICATION  MATRIX; 

61. Pl_ 

G1  24  1 

PI  0 25 

G2  24  1 

P2  10  15 


F-LEVEL  FOR  DELE  TION: 


0.005 


INITIAL 

FINAL 

U-STAT 

0.43226 

Approximate  -F 
6a.  04376 

0.19181 
— *0. 1 9640 

Red 

Green 

Blue 

Yellow 

White 

Purple 

Blue-G 
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THE  BCI  FACILITY 


Following  is  a brief  description  of  the  Brain  Computer  Interface  Labora- 
tory as  implemented  during  the  period.  An  addition  of  a 16K  XDS  930,  a 
random  access  drum,  printer,  two  7-track  magnetic  tape  drives  and  a card 
reader  to  the  BCI  experiment  loop  was  being  made  at  the  end  of  year  1974, 
and  details  on  the  new  hardware  and  its  role  will  appear  in  subscnuent 
Reports. 

OPERATIONAL  SYSTEM  DESCRIPTION 


Three  computers  are  involved  in  the  Brain  Computer  Interface  Project,  the 
IMLAC  PDS-1,  the  XDS  920  and  the  IBM  360/91.  High  speed  parallel  data 
links  connect  the  three  computers  through  a common  interface  box.  The 
interface  box  allows  the  three  computers  to  talk  to  each  other  by  using  a 
complex  heirarchy  of  interrupt  control. 

The  subject  is  placed  in  a sound  proof  room  (see  Figs  29  and  30  ).  The 
room  also  shields  the  experiment  from  magnetic  and  electric  fields  in- 
cluding radio  frequencies  as  well  as  allows  referenct  to  a single  point 
to  ground.  The  experimenter  sits  in  the  adjoining  room  in  front  of  the 
contro  terminal.  The  EEG  signals  are  preamplified  and  monitored  and 
plotted  using  a standard  EEG  recording  machine.  The  amplified  signals  are 
also  sent  to  a 50  channel  analog  to  digital  converter  which  converts  and 
feeds  the  signals  to  the  XDS  920  computer.  This  XDS  920  acts  as  self 
driven  data  input  controller  and  real  time  scheduler  for  the  experiment. 

It  also  performs  real  time  processing  functions  used  in  the  creation  of 
the  EEG.  These  data  sets  contain  complete  information  regarding  the 
experiment. 

The  other  two  computers  are  the  IBM  360/91  and  the  IMLAC  PDS-1.  The 
IBM  360/91  provides  the  main  computing  power  for  the  experiments.  Most 
of  data  reduction  and  data  analysis  is  performed  on  the  IBM  360/91  due 
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to  the  availability  of  BMD  statistical  analysis  software,  its  exceptionally 
large  core  memory  and  its  excellent  number  crunching  abilities. 

The  laboratory  has  a direct  access  post  into  the  360  through  which  one 
can  read  from  or  write  into  a dedicated  120K  of  core  memory. 

A MONITOR  program  which  controls  the  data  flow  as  well  as  the  processing 
schedules  has  been  implemented  on  the  360/91.  The  MONITOR  is  initiated 
by  using  the  campus  time  shared  URSA  system  through  a CCI  CRT  terminal 
installed  in  the  laboratory. 

Among  the  various  analysis  softward  on  the  360  are  pattern  recognition 
programs  to  study  and  quantify  the  discriminabil ity  of  data  belonging  to 
varying  stimuli. 

The  IMLAC  PDS-1  is  a graphic  computer  which  drives  a slave  display  located 
in  a window  in  front  of  the  subject  to  provide  visual  feedback  display. 

The  IMLAC  PDS-1  is  primarily  used  in  pattern  related  evoked  response  exper- 
iments, including  variation  of  horizontal  and  vertical  lines,  faces,  mazes 
etc. 

DATA  ACQUISITION 
EEG  RECORDING  MACHINE 


The  EEG  recording  machine  is  a Grass  Medical  Instrument  Model  6.  The 
instrument  uses  plug  in  type  units  with  color  coded  pushbuttons  for 
electrode  selection  and  regulated  voltages,  hybrid  amplifiers  (with  60 
cycle  filters). 

The  electrode  board  with  the  International  Numeric  System  has  23  numbered 
jacks  for  all  possible  positions  of  electrodes. 

A master  gain  and  reduction  control  switch  can  adjust  all  channel 
amplifications. 
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ELECTRODES 


The  electrodes  used  in  the  experiments  are  surface  disc  silver  chloride 
electrodes.  A compound  electrolyte  is  used  to  provide  contact  with  the 
scalp  and  to  secure  the  electrode.  The  area  of  the  scalp  is  first  pre- 
pared with  alcohol  and  abrased  slightly  to  eliminate  skin  oils  and  de- 
crease scalp  resistance. 

ANALOG  INPUT  TO  DIGITAL  DATA  CONVERSION 

A/D  CONVERTER 

A SDS  Model  AD  20-11  made  up  of  silicon  semiconductor  discretes  which 
is  capable  of  30000  12-bit  conversions  per  second  is  used,  ine  converter 
is  accurate  up  to  99.95%,  accepts  up  to  10  volts,  and  has  a Msual  dis- 
play of  the  latest  conversion.  An  aperture  of  up  to  100  msec,  can  be 
realised  with  an  optional  sample  and  hold.  The  average  conversion  time 
is  33.3  microseconds.  The  converter  has  a self-contained  reference 
voltage  generator. 

MULTIPLEXERS 


One  16  channel  and  one  32  channel  basic  SDS  MU31  multiplexers  are  installed. 
D/A  CONVERTER 

Model  DA  30-15  is  used  to  convert  12  bits  of  parallel  digital  data. 

The  converter  is  used  to  drive  the  analog  plotter  and  a CRT. 

IMPLEMENTATION  OF  50  CHANNEL  A/D  CONVERSION 

The  analog  input  is  converted  to  digital  data  using  a complex  of  4 A/D 
converters,  2 multiplexers  and  3 D/A  converters.  The  analog  input  data 
is  sent  to  two  of  the  ADC  directly  whereas  the  other  two  ADE's  work  with 
the  two  multiplexers  by  processing  two  channels  at  a time.  The  converted 
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data  from  ADC  #1  is  sent  to  the  left  half  of  the  24  bit  word  and  the  con- 
verted data  from  ADC  #2  is  sent  to  its  right  half.  The  ADC  outputs  in  a 
(10  data  bits  + 1 sign  bit)  register.  Thirty-two  channels  would  be  pro- 
cessed in  a little  over  960  microseconds.  A sampling  time  is  usually 
one  to  eight  milliseconds. 

THE  IMLAC  PDS-1 

The  IMLAC  PDS-1  is  a graphical  terminal  and  minicomputer  which  gen- 
dtps  pattern  display  stimuli  for  the  subject.  It  is  a 8K  dual  pro- 
cessor in  which  the  two  processors  share  the  same  core  memory.  One 

of  the  processors  is  a conventional  16  bit  minicomputer,  whereas  the 
other  is  a display  processor. 

The  data  processor  has  a cycle  lime  of  2 ^sec.,  J1  ect  addressing  of 
2K  and  indirect  addressing  up  to  32K.  It  has  one  general  purpose  regis- 
ter plus  one  LINK  bit.  It  allows  immediate  instructions  and  has  parallel 
16  bit  I/O  up  to  a maximum  of  100K  words  per  second  with  one  level 
priority  Interrupt.  It  allows  power/fail  restart  for  core  as  well  as  for 
display  information.  Besides  the  usual  set  of  instructions  it  has  instruc- 
tions for  turning  the  display  processor  on,  enable  and  disable  interrupts, 
high  voltage  synchronization,  continuous  2 msec  intensification,  TTY, 
paper  tape  and  light  pen. 

THE  IMLAC  DISPLAY  PROCESSOR 

The  display  processor  oses  its  own  instruction  set  and  is  capable  of 
decoding  Instructions  as  well.  It  steals  memory  cycles  from  the  mini- 
computer when  necessary.  It  can  function  in  two  modes;  the  Increment  mode 
and  the  Processor  mode.  The  later  allows  it  to  decode  a small  instruction 
set,  whereas  the  former  allows  it  to  generate  screen  displays  by  decoding 
memory  words  into  beam  movements  and  corresponding  Intensification. 

The  long  vector  hardware  allows  it  to  draw  lines  up  to  full  length'  of 
the  screen  as  well  as  dotted  lines k 


151 


) 


The  display  in  the  sound  proof  subject's  room  operates  in  a slave  mode 
to  the  display  processor  providing  pattern  display  stimuli.  The  IMLAC  PDS-1 
generates  displays  and  puts  them  on  the  screen  before  the  experimenter  on 
the  master  display  Independently  of  the  slave  display.  The  software  for 
IMLAC  Includes  a graphic  design  language  (GRAL),  SPACEWAR  program  and  face- 
qeneratlng  programs  which  are  scheduled  to  be  u-,ed  in  experiments. 

THE  XDS  920 

The  function  of  this  computer  is  to  acquire  data  from  the  analog  to  digital 
converter  and  to  act  as  a scheduler  of  the  real  time  experiment  control. 

The  DHS  (Data  Handling  Supervisor)  controls  the  real  time  collection  of 
data,  its  analysis  and  display.  It  converses  with  the  other  computers 
by  means  of  a 24  bit  parallel  data  link  which  can  operate  at  the  rate  of 
40K  bit  per  second  through  an  Interface  box. 

It  a computer  sultaole  for  on-line  and  real  time  systems  that  require 
high  processing  rates,  complex  I/O,  and  control  functions.  However  It 
has  a limited  memory  of  4096  word;.  The  computer  can  handle  up  to  65K 
words  per  second  In  data  transmission  operation.  It  Is  a binary  computer 
with  a ? s complement  facility  and  has  a memory  cycle  time  of  8 msec.  It 
allows  Interrupts  and  uses  a real  tine  clock  which  can  be  programmed  to 
produce  timing  pulses  from  500  cps  to  8000  cps  to  provide  interupts. 

The  computer  uses  a XDS  MAGPAK  unit  which  facilitates  programming  and 
operation.  It  can  be  used  to  store  library  compilers  and  assemblers 
and  could  be  ooo*strapped  in  with  ease.  MAGPAK  has  two  tape  cartridges, 
each  tape  having  a capacity  of  2.8  million  bytes  of  6 b'ts  plus  parity  each. 
The  transfer  rate  Is  1500  characters  per  second.  It  is  used  to  hold: 

! ibrary.  Temporary  Storage,  Object  programs  and  Source  programs.  A XDS 
9150  card  reader  Is  also  used  which  reads  cards  up  to  100  card-  oer  minute. 
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THE  IBM  360/i)l  I I 

This  campus  machine  provides  the  main  computing  power  due  to  Its  excep- 
tional large  core  memory  and  fast  number  crunching  abilities  The 
computer  Is  used  through  a CCI  CRT  terminal  located  in  the  laboratory 
under  the  campus-wide  time-shared  URSA  system. 

This  computer  houses  the  main  data  reduction  and  analysis  software  which 

can  be  used  through  the  MONITOR  program  which  is  also  installed  on  the 

360  and  controls  the  data  flow  between  the  lab  and  the  360  as  well  as 
the  data  analysis. 
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COMPUTER  INTERFACE  FACILITY 


PREAMPLIFIERS 
& FILTERS 


A/D 

ANALOG  L 
TO  DIGITAL 
CONVERSION 


DIGITAL 

COMPUTER 


DIGITAL 
TO  ANALOG 
CONVERSION 


SYSTEM  DIAGRAM  OF  A BRAIN  COMPUTER-INTERFACE 
FIG.  32 


i 


55 


